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ABarRACT 


Ballistic  resistance  of  thin  metallic  plates,  both  single  plates  and 
laminates,  to  projectile  Impact  Is  investigated.  Two  mathematical  models  are 
formulated,  each  emphasizing  a  critical  aspect  of  the  deformatlonal  process. 

The  first  model  considers  failure  by  high  tensile  or  shear  stress  in  the 
transient  stress  wave  that  propagates  across  the  plate,  while  the  second  model 
considers  failvire  by  sustained  plastic  shear  flow. 

Engineering  models  of  fracture  and  flow  under  high  speed  impact  are 
evaluated,  using  the  continuum  theory  of  dislocations.  These  models  are 
used  In  the  analysis  of  the  modes  of  failure.  Critical  inpact  velocities 
for  failure  are  calc\ilated  and  characterized  as  a  function  of  dynamic  material 
properties. 

The  resvilts  of  this  study  Indicate  that  there  Is  a  need  for  research  on 
the  determination  of  high-speed  properties  of  materials.  Improved  armor 
materials  can  be  realized  throxogh  materials  research  based  upon  theories 
developed  herein. 

Further  analytic  research  on  stress-wave  propagation  and  failtire  in 
materials  other  than  metals  is  advised.  Both  the  present  analytic  work  and 
future  studies  should  be  used  to  formulate  and  carry  out  an  experimental 
program. 
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Section  1 


IHTIfflrUCTION 

Throvighout  the  annals  of  ancient  history  there  are  many  references  to  the 
arms  and  armor  which  great  armies  bore  in  battle.  Horner  makes  reference  to 
arms  and  armor  made  of  bronze  and  occasionsJ-ly  of  iron.  Body  armor  at  first 
consisted  merely  of  shields  and,  a  little  later,  helmets.  It  was  not  too  long 
however  before  breast  plates  and  shin  guards  were  in  vogue.  During  the  middle 
ages,  complete  body  armor  was  the  standard  uniform  of  the  equestrian  knight; 
indeed,  even  his  horse  was  protected  by  form-fitting  metal  plates.  It  did  not 
take  too  long  a  period  of  time  before  the  suits  of  armor  became  vehicles  for 
artistic  expression  in  addition  to  their  protective  purposes.  Exan^les  of  this 
form  of  art  are  now  the  prizes  of  many  art  collectors.  Examination  of  these 
relics  of  a  bygone  era  will  show  why  they  were  reserved  only  to  horse-borne 
warriors.  The  wearer  was  not  very  mobile.  It  was  long  after  the  use  of  such 
armor  became  general  that  it  was  doomed  to  become  obsolete.  This  was  brought 
about  first  by  the  efficiency  of  the  English  long  bow  and  finally  by  the  coming 
of  age  of  firearms. 

As  firearms  came  into  Increasing  use  and  brought  about  a  new  concept  of 
warfare,  one  which  depended  on  a  mobile  infantryman,  body  armor  became  increasingly 
unpopvilar.  The  greater  and  greater  muzzle  velocities  which  firearms  were  becoming 
capable  of  were  decreasing  the  feasibility  of  individualized  protection.  In 
order  to  stop  a  high-velocity  bullet,  the  armor  woxald  have  had  to  be  so  thick 
that  it  woiild  be  impossible  for  a  man  to  carry  and  still  be  mobile. 

The  concept  of  body  armor  has  not  been  abandoned,  however.  The  quest  for  an 
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lndlvld\ial  protective  shield  that  will  he  effective  against  modem  arms  and  yet 
leave  its  wearer  coo5>letely  mobile  still  continues.  Up  to  very  recent  tlmes> 
this  quest  has  been  \inldirectlonal  In  approach  emphasizing,  almost  exclusively, 
the  experimental  viewpoint.  It  does  not  come  as  a  conqplete  surprise  that  the 
results  of  this  approach  have  not  been  overly  gratifying. 


Throughout  the  history  of  science  there  have  been  many  proponents  of  the 
strictly  experimental  approach  as  well  as  the  strictly  contemplative  or  theoretical 
approach  to  given  problems.  The  ultimate  lesson  to  be  learned  from  this  is  that 
neither  is  completely  adeqiiate  by  itself.  And  so  it  follows  that  a  wedding  of 
these  two  approaches  is  the  best  way  to  solve  a  given  practical  problem. 


Achieving  such  a  iinlon  is  not  an  easy  task, however.  There  is  generally 
a  wide  gap  existing  between  theory  and  experiment.  This  is  mainly  due  to  the 
difference  of  interests  of  the  theoretician  and  experimenter  in  addition  to 
at  least  a  modlc^xtn  of  mutual  distrust.  The  present  study  is  directed  toward 
the  development  of  an  adequate  theory  that  will  enable  the  experience  of  existing 
experimental  studies  to  be  given  a  more  meaningful  interpretation.  In  addition, 
and  perhaps  more  important,  the  goeGL  of  the  present  theoretical  endeavor  is  to 
provide  a  more  fruitfiil  basis  for  future  experimentation. 


Thus  the  present  study  is  intended  to  be  the  first  theoretical  link  in  a 
balanced  overall  program  directed  toward  the  design  of  personnel  armor  that  will 
be  effective  against  modern  firearms. 


1.1  Objectives 

In  view  of  the  state  of  the  art  of  body-armor  design  that  was  cvirrent  when 
the  present  study  began,  it  was  not  difficult  to  determine  the  objectives  of  the 
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program.  Indeed,  the  ohjectlves  were  dictated  by  the  circumstances^) 

First  of  all,  a  theory  must  be  developed  to  describe  what  occurs  when  a 
bullet  strikes  a  piece  of  armor.  This  theory  must  be  based  on  the  basic  concepts 
of  the  physical,  laws  which  govern  the  behavior  of  materials  under  the  conditions 
of  impact.  The  role  that  the  properties  of  the  materials  involved  play  must  be 
determined.  The  study  must  ascertain  which  properties  or  combinations  of  properties 
have  the  greatest  bearing  on  the  adequacy  of  the  armor  to  do  its  job.  The 
theory  developed  must,  wherever  possible,  take  into  account  the  existing  res\ilts 
of  the  experimenters  in  order  that  there  be  communication  with  the  realities 
of  experience. 

Finally,  the  theory  must  point  out  the  path  for  future  experiments.  It 
must  answer  the  question  of  the  experimenter,  "What  area  have  I  failed  to  explore?" 
Then,  when  the  experimenter  has  responded  with  new  experiments  and  the  theoretician 
has  joined  him  in  reaching  a  compatible  answer  to  the  problem  of  penetration,  the 
designer  may,  with  knowledge,  approach  his  task. 

1.2  Previous  Work 

As  may  be  noted  from  the  title  of  this  report,  there  has  been  a  previous 
volume  dealing  with  the  mechanics  of  penetration.  It  was  the  task  of  Voliome  I 
to  break  the  ground  for  the  theory  of  penetration.  Outside  of  a  few  voices 
that  had  cried  in  the  wilderness  (notably  Bjork^^\  and  Hopkins  and  Kolsky^^"^) 
there  had  been  little  done  in  the  way  of  theoretical  work  on  the  problem. 

Volume  I  gave  the  res\alts  of  the  first  year's  study  of  the  mechanics  of 
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penetration.  While  this  preparatory  work  took  cognizance  of  the  aspects  of 
plastic  flow  involved  in  the  problem,  it  did  not  attempt  to  solve  the  problem 
in  terms  of  plasticity.  This  woilLd  have  been  too  great  a  bite  for  an  area  in 
which  so  little  work  had  been  done.  Instead,  the  phenomena  were  investigated 
from  the  viewpoint  of  elasticity  theory,  using  this  approach,  a  framework 
for  more  advanced  study  was  constructed.  Qualitative  descriptions  and  quantitative 
boundaries  were  established.  The  basic  mechanism  of  the  phenomena  could  be 
worked  out  without  becoming  entangled  in  the  intricate  niceties  of  plastic  theory. 

This  voliome,  based  on  the  ground  work  of  Volume  I,  now  attacks  the  f\all- 
blown  problem  and  considers  the  effects  of  plastic  flow.  It  investigates  the 
effects  of  plastic  deformation,  fracture,  and  the  two  types  of  failure  that  may 
occur,  tensile  failure  or  shear  failure.  It  considers  the  combined  effect  of 
these  two  types  of  failure  when  applied  to  a  given  problem.  The  theory 
developed  will  be  applied  to  the  penetration  of  laminated  metallic  plates. 
Recommendations  for  future  work  in  the  area  of  armor  penetration  are  made 
based  on  the  conclusions  that  may  be  drawn  from  this  study. 

1.3  The  Physical  Phenomena  of  Fracture 

When  a  target  is  Impacted  by  a  projectile,  a  number  of  physical  phenomena 
may  occur,  depending  on  the  magnitude  and  direction  of  the  impact  velocity, 
and  on  the  materials  and  geometries  of  the  target  and  projectile.  For 
prescribed  impact  conditions,  the  relative  importance  of  these  phenomena  will 
determine  the  nature  of  the  impact  and,  therefore,  such  qroantitative  aspects 
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as  momentvim  transfer,  partition  of  the  projectiles,  kinetic  energy  into  other 
forms,  and  the  type  and  extent  of  fracture.  For  a  specific  target  and  projectile, 
the  impact  velocity  alone  will  determine  the  character  of  the  Impact,  and, 
consequently,  it  is  possible  to  divide  the  general  impact  problem  into  a 
number  of  regimes  that  depend  on  the  impact  velocity.  In  each  of  these  regimes, 
a  different  phenomenon  or  combination  of  phenomena  is  important  in  controlling 
the  physical  processes  that  occur.  There  is,  of  course,  considerable  overlap 
between  the  regimes,  but,  nevertheless,  they  are  a  valuable  way  of  describing 
the  various  kinds  of  impact.  In  order  of  increasing  impact  velocity,  the 
regimes  are : 

1.  Elastic  Impact 

2.  Plastic  Impact 

3.  ^drodynamic  Impact 

4.  Superset smic  Impact- 

5 .  Explosive  Impact 

The  first  of  these  regimes,  elastic  impact,  is  characterized  by  the  fact 
that  the  stresses  produced  by  impact  in  the  target  and  in  the  projectile  are 
not  above  the  elastic  limit,  and  the  impact  process  is,  consequently,  fully 
reversible.  For  this  type  of  impact,  the  kinetic  energy  of  the  impinging 
projectile  is  usually  converted  almost  entirely  into  elastic  strain  energy 
in  the  target  and  projectile.  Elastic  stress  waves  will  be  generated  in  the 
target,  consisting  of  primary  (or  dilatational)  waves,  secondary  (or  distortions! ) 
waves,  and  Rayleigh  surface  waves.  However,  the  duration  of  impact  at  the 
low  velocities  associated  with  elastic  impact  is  so  long  that  there  is  time 
for  many  reflections  of  the  waves  within  the  target  and  projectile,  and  a 
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quasi- static  stress  distribution  is  effectively  established^  with  very  little 
of  the  total  energy  converted  into  the  kinetic  energy  of  elastic  waves.  A 
reasonably  acciorate  description  of  this  type  of  Impact  can  be  obtained  by 
presuming  that  the  entire  kinetic  energy  of  the  incoming  projectile  is  con¬ 
verted  into  the  strain  energy  of  a  static  stress  distribution  as  given,  for 
instance,  by  the  Hertz  formulation  for  contact  stresses. 

When  the  impact  velocity  is  increased  beyond  the  critical  velocity 
where  the  elastic  limit  of  either  the  target  or  projectile  is  exceeded,  the 
stress  distribution  will  be  partly  elastic  and  partly  plastic,  which  is  typical 
of  the  regime  of  plastic  impact.  The  critical  velocity  that  separates  the 
elastic  and  plastic  regimes  depends  markedly  on  the  materials  and  also  on  the 
geometries  of  the  target  and  projectile.  It  has  been  found,  for  example, 
that  the  critical  velocity  for  a  l/4-inch  diameter  steel  ball  impacting 
normally  onto  a  large  steel  plate  is  only  about  0.15  m/sec^^^\  For  the  normal 
Impact  of  two  flat-ended  steel  cylinders,  however,  the  critical  velocity  is 
in  excess  of  6o  m/sec.  The  criticsG.  velocity  for  plastic  impact  also  depends, 
of  course,  on  the  elastic  limit  of  the  materials  of  the  target  and  projectile, 
and  if  the  materials  behave  in  brittle  rather  than  a  ductile  manner,  no 
plastic  regime  will  exist. 

The  regime  of  hydrodynamic  impact  is  characterized  by  stresses  that  are 
far  in  excess  of  the  yield  strength  of  the  materials,  with  the  resiolt  that, 
for  a  part  of  the  time  of  Impact,  the  materials  flow  freely.  The  sequence 
of  such  an  impact  is  depicted  schematically  in  Figure  1.1^^^\  Immediately 
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Figure  1.1  -  Sequence  of  Qy'drodynamlc  Tmjftryh 
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after  contact,  a  region  of  extremely  high  conqjressive  stress  is  created  near 
the  projectile-target  Interface,  and  a  con5>resBlve  wave  is  initiated.  Shortly 
thereafter,  the  material  of  hoth  target  and  projectile  begin  to  flow,  and  the 
compressive  wave  and  shear  wave  propagate  into  the  body  of  the  target.  The 
flow  of  the  projectile  material  continues  until  the  projectile  is  entirely 
dissipated,  and  the  flow  in  the  target  continues  for  a  somewhat  longer  period 
because  of  the  very  high  flow  velocities  induced  by  the  projectile.  Finally, 
the  crater  formed  by  the  impact  contracts  slightly  because  of  elastic  (and 
some  plastic)  rebound.  This  type  of  impact  can  be  described  theoretically 
by  hydrodynamic  calcTilatlons,  which  assume  that  the  materisG-s  behave  as  fluids 
with  no  effective  resistance  to  distortion.  Such  calculations  are  accurate 
in  determining  the  flow  patterns  and  stress  distribution  when  the  stress  levels 
are  well  above  the  yield  strength  of  the  materials,  but  they  are  not  accurate 
in  predicting  the  later  stages  of  the  intact  (afterflow  and  rebound)  where 
the  yield  strengths  are  important. 

The  sequence  of  events  that  comprise  an  impact  in  the  superseismic 
regime  is  shown  in  Figure  1.2.  The  Initial  stage  of  high  compressive  stress 
(in  the  megabar  range)  is  similar  to  the  corresponding  state  of  hydrodynamic 
impact,  except  for  considerably  higher  stress  levels.  The  next  stage,  however, 
is  different  in  that  no  compressive  wave  will  propagate  away  from  the  expanding 
cavity  into  the  body  of  the  target  material,  because  the  cavity  is  expanding 
at  superseismic  speed.  The  lateral  and  rearward  flow  of  the  materials  continues 
in  the  same  general  manner  as  for  hydrodynamic  impact,  except  for  the  higher 
flow  velocities.  Most  inqjortantly,  the  flow  can  occur  only  within  the  super- 
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Figure  1.2  -  Sequence  of  SuperselBinlc  Inqpact 
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selsmlcally  expanding  cavity,  and  since  the  dcfwnward  penetration  of  the  pro¬ 
jectile  is  limited  hy  the  seismic  velocity,  the  flow  of  material  is  general. 
thro\ighout  a  hemispherical  region.  As  a  result,  the  final  crater  of  a  suijer- 
seismic  impact  is  typically  hemispherical. 

The  regime  of  explosive  Impact  is  characterized  by  impact  velocities 
high  enoiigh  to  vaporize  material  of  the  target  and  projectile.  Under  these 
conditions,  the  Impact  problem  Is  akin  to  a  surface  explosion,  wherein  the 
momentum  transfer  and  stresses  are  produced  by  the  high  pressure  of  vaporized 
materials. 
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1.4  Scope 

In  Section  2^  a  qualitative  description  of  plastic  deformation  Is  presented. 
The  mathematical  form  of  the  constitutive  equation  is  given  and  the  evaloiation 
of  specific  models  Is  made.  An  Interpretation  of  dynamic  yield  strength  Is 
given  and  discussed.  The  material  contained  in  this  section  represents 
original  work  in  this  field. 

The  subject  of  fracture  is  treated  in  Section  3.  The  Griffith  theory 
of  fractiire  is  explained  and  the  growth  of  macro-holes  is  discussed  for  the 
special  cases  of  interest.  The  relationship  of  tlme-to -fracture  to  the 
problem  of  penetration  is  then  explored. 

Sections  4 #5’  &  ^  deal  with  the  subject  of  failvire  of  both  single  and 
laminated  plates.  Both  tensile  and  shear  failure  and  a  combination  of 
these  are  investigated.  The  equations  of  motion  and  their  solutions  are 
presented  and  numerical  results  arfe  given.  Conclusions  and  recommendations 
are  presented  in  Section  7- 
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Section  2 


ELASTIC  DEFORMATION  OF  MCTAL  UNDER  TRANSIENT  LOADS 

In  this  prohlem  of  impact  deformation  and  penetration  of  metallic  plates 
by  projectiles,  the  stresses  in  the  plate  often  exceed  the  static  proportional, 
limit.  When  the  intact  stress  exceeds  this  level  the  constitutive  relations 
between  the  stress  and  the  strains  are  markedly  altered. 

The  basic  mechanism  for  plastic  deformation  in  metals  is  by  the  movement 
and  formation  of  dislocations  under  stress.  The  physical  backgrovind  for  this 

(i4^ 

mechanism  was  discussed  in  Vol.  I 


2.1  Constitutive  Equation 

The  form  of  the  constitutive  equations  for  a  plastic  flow  in  crystalline 

and  ploycrystalline  material  specimens  has  been  investigated  quite  thoroughly 

(l2)(25) 

in  view  of  microscopic  mechanism  of  plastic  flow'  ^  , 


In  an  inelastic  material,  the  rate  of  total  strain  may  be  decomposed 
into  two  parts,  the  strain  rate  corresponding  to  the  elastic  distortion 
or  strain  of  the  sample  and  to  the  plastic  strain  rate  or  rate  of  permanent 
deformation. 


t 


IJ 


(2.1) 


where 


£  =  total  strain  = 

J 


.  .  du  ^ 

4  (—  +  — ^) 


£  = 
IJ 

e  P  = 
ij 

(•)  = 


elastic  strain 

plastic  or  permanent  strain 

— (  ) 
dt''  ^ 
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The  elastic  strain  is  related  directly  to  the  stress.  If  the  linear  elastic 
relationship  (Hooke's  Law)  holds,  this  component  of  strain  rate  is  proportional 
to  the  stress  rate. 


.  E 
'ij 


(2.2) 


°ijk£~  elastic  constants 

If  the  mechanism  of  plastic  deformation  Is  the  passage  of  a  n\imber  of 

dislocations  through  the  sample  In  response  to  the  applied  stress,  the  expression 

(12) 


for  the  plastic  strain  rate  component  may  be  written 
t 


where 


'IJ  “'rq 


(2.3) 


'rqj 


a 


iq 


V 


1  rqj  even  permutation  of  1,2,3 
-1  rqj  odd  permutation  of  1,2,3 
0  otherwise 

is  the  number  of  dislocations  per  unit  area  having  its  Burger's 
vector  in  the  1-direction  and  having  its  axis  in  the  q-directlon 

is  the  velocity  component  of  the  dislocation  in  the  r  direction. 


The  siimmation  convention  over  repeated  indices  is  used.  The  bar  over  the 
whole  expression  indicates  that  the  average  is  taken  over  all  possible 
orientations  of  the  dislocation  in  the  plane  normal  to  its  axis. 


In  this  expression,  the  velocity  of  the  dislocation  is  a  function  of  the 

shear  stress  in  the  plane  normal  to  the  axis  of  the  dislocation.  The  increase 

(12) 

of  velocity  is  not  linear  function  of  stress.  Fugelso  has  worked  out 
a  theoretical  expression  for  this  velocity  in  terms  of  this  shear  stress. 

_  M 

V  =  Vjj  e  slnh7T  1^1*^  1 

V  =  Vjj  It  |>  1 
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where 


AH 

T  Is  the  solution  to  slnh^T  = 
c 

A1  =  activation  energy  for  self-diffusion 
R  =  gas  constant 
T  =  absolute  ten^jerattire 

The  ntamber  of  dislocations  Is  incorporated  In  the  definition  of 
This  n\jmber  may  be  constant  or  may  increase  by  creation  or  regeneration 
mechanisms. 

The  behavior  of  several  singer  problems  will  now  be  investigated. 

These  cases  will  be  chosen  for  their  diagnostic  value,  first  In  checking  the 
validity  of  the  model  for  plastic  flow  for  form,  second  to  determine  basic 
methods  for  determining  the  niamerical  value  of  the  constants  thro\;igh  controlled 
experiments. 

Consider  the  case  of  uniform  deformation  of  a  rectangular  block  of 
material • 


The  remaining  constitutive  eq\iations  may  be  written 


t  =1 

XX  E 

»xx 

-  |(»yy 

+ 

*zz) 

+ 

t  p 

XX 

t  -  -p 

yy  E 

*yy 

■*  -E^^xx 

•f 

»zz) 

+ 

fe  P 
yy 

^zz  =  4 

»zz 

-  -i^^xx 

+ 

+ 

t  P 
zz 

(2.4) 

with 

E  =  Young's  Modulus 
V  z  Poisson's  Ratio 
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since  the  stemdard  tensile  and  cccnpresslve  test  apparatus  set  up  a 
state  of  biaxial  stress  —  this  model  will  be  restricted  to  biaxial  symmetry 


yy  2z 

&  =  & 
yy  zz 


The  general  form  of  plastic  con^xDnents  of  strain  have  been  evaluated 
for  the  case  of  biaxial  stress' 


(2.5) 


£  P  =*  £  P  =  VjMh 
YY  zz  D 


(2.6) 


where 


AH 

o 


Vjj  =  Vjj  e  slnh7.T 


V  =  V 
''d  '^D 


'a  2  +  a 

yy  zz  > 


[r|  < 


T  >  \T, 


is  given  by' the '  solution  of  slnh7T' 


M/rt 


Vp  is  the  dislocation  velocity 


Vjj  is  the  maximum  dislocation  velocity 

N  is  the  number  of  dislocations  per  unit  area 
b  is  the  Burger's  vector  of  the  dislocation 


In  the  most  elementary  testing  machine  the  lateral  stresses  are  zero 


‘^yy  "^zz  “  ° 


(2.7) 
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II 


Thus 


"2  ^xx  only  the  first  equation  is  left 


t 


=  2Vr^bN 


(2.8) 


Inserting  the  value  of  the  dislocation  velocity,  the  constitutive  equation  is 


'xx  "  E 


-  2V  °  e"^  Blnh„ 


a  <  2t 
XX  c 


^xx  =  ^ 


a  >•  2t 
XX  c 


(2.9) 

Rather  than  evaluating  this  pair  of  equations  with  the  discontinuous  form, 
we  shall  replace  the  expression  for  the  dislocation  velocity  with  an  analytic 
expression  over  the  range  0  ugl  which  is  very  close  approximation 


=  V 


Vj,  =  -  Vp°  e-AA' 


O  <  r  < 


<  T  <  0 


The  comparison  of  the  two  expressions  is  given  in  fig.  2.1  .  Only  in 

the  neighborhood  of  zero  shear  stress  is  there  any  marked  deviation  and  in 
this  region  the  contribution  from  this  term  is  negligible  there  anyway. 


The  coefficient  A  absorbs  the  temperature  behavior 
*  ET 


(2.10) 
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Dislocation  Velocity^!-'  fps 


\ - 1 - 1 - 1 - 1 - 1 - 1 - 1  r 

0  12 


Shear  Strain  (kbar) 


1 - r 
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This  form  is  preferable  since  ve  will  be  dealing  only  vlth  solutions 
of  problems  whose  stresses  have  the  same  sign  during  the  process  of  deformation 


Four  cases  will  now  be  considered 

a)  Constant  strain  rate 

b)  Constant  strain 


Constant  nusiber  of  dislocations 


c)  Constant  stress  j 

d)  Consitant  strain  rate  -  Increasing  nuiriber  of  dislocations 
We  may  drop  the  subscript  xx  as  no  confusion  will  arise. 

Case  a)  Constant  strain  rate 


This  problem  will  be  solved  for  a  constant  strain  rate 
£  =  C  =  constant 


Inserting  this  Into  the  constitutive  equation 


(2.11) 


This  equation  Is  solved  n\amerlcally  by  the  Runge-Kutta  method  for  several 
values  of  C,  the  applied  strain  rate.  The  values  of  C  range  over  several 
orders  of  magnitudes  from  10”Vsec  to  10”^/sec.  (Fig.  2.2) 

It  Is  found  that  the  stress  Increases  linearly  with  strain  at  the 
beginning  and  then  very  rapidly  changes  slope.  This  sharp  bend  occurs  at  a 
stress  which  is  far  below  the  maximum  possible  plastic  flow  rate.  If  Its 
applied  strain  rate  Is  less  than  the  maximum  plastic  strain  rate/  the  stress 
approaches  a  limiting  value  which  Is  given  by 

PA 

=  Iri  (2.12) 
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The  plastic  flow  Is  able  to  absorb  the  entire  applied  strain.  This  value 
of  the  stress  may  be  considered  to  be  the  dynamic  yield  strength  of  the 
material  at  that  strain  rate.  If^  on  the  other  hand>  the  critical  strain  rate 
exceeds  the  maximum  flow  rate^  the  plastic  flov  cannot  respond  fast  enough 
to  account  for  all  the  strain  and  the  stress  must  Increase  Indefinitely. 


Case  b)  Constant  strain^  suddenly  applied  and  held 


If  Is  the  value  of  the  strain  at  t  =  0  and  Is  the  stress  at  t  =  0 
o  o 

=  Ee^  since  no  flow  has  occured 
o  o 

and 

t  =  0 


Thus 


da 

■3^ 


-2A/a 

e 


0 


(2.13) 


Subject  to  the  InltleO.  conditions 

a  =  a  at  t  =  0 
o 

the  stress-time  history  is  plotted  in  Figure  2.3, 

It  is  seen  that  the  stress  decays  from  its  initial  value  very  rapidly  if  the 
initial  stress  is  high  and  very  gently  if  the  Inltied.  stress  is  small. 


Case  c)  Constant  stress 

At  time  t  =  0,  the  stress  Ib  suddenly  applied  and  maintained 


a  =  const  =  a. 


thus  the  constitutive  equation  is 

_  -2k/ a 


(2.11^) 
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The  solution  to  this  equation  Is 


(2.15) 


The  strain  starts  at  Its  elastic  value  and  Increases  linearly  with  time. 
The  rate  of  Increase  of  strain  approaches  a  limiting  value,  as  the  stress 
becomes  very  large. 


Case  d)  Constant  strain  rate  -  Increasing  number  of  dislocations 

In  materials  the  munber  of  dislocations  may  increase  by  a  variety  of 

(13  ) 

microscopic  mechanisms,  Fugelso  has  evaluated  thf  ■‘havior  of  the  stress- 
strain  curves  if  the  number  of  dislocations  increases  by  the  Prank-Read  dis¬ 
location  log  mechanism.  In  terms  of  macroscopic  variables  the  constant 


term  in  the  expression  for  strain  rate  is  replaced  by 


(e,i6) 


The  equation  for  stress- strain  at  constant  rate  becomes 


C  =  —  +  ^  S  e 

dt  ^  o 


(2.17) 


This  equation  was  solved  nximerlcally  for  several  velues  of  the  applied  strain 
rate,  fig.  2.I4.. 


In  this  case  the  stress- strain  curve  is  altered.  The  stress  rises  rapidly 
along  its  elastic  veilues,  passes  through  a  maximum  as  the  plastic  flow  becomes 
dominant,  decays  rapidly  and  finally  reaches  a  plateau  where  the  decay 
is  very  slow. 
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Strain  rate^  sec 


Figure  2,5  -  Eiyuamlc  Yield  Strength  as  a  Function  of  Strain- 
Bate  for  a  T1  2.5  Sn  5  A1  Target 
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2.2  Hynamlc  Yield  StreBses 


Consider  the  resulting  stress-strain  curves  at  constant  strain  rate. 

For  any  applied  strain  rate  below  the  naximum  rate,  the  stress  rapidly  levels 
off  to  a  constant  value  whlcn  is  characteristic  of  the  strain  rate.  This 
stress  is  then  maintained  throughout  the  remainder  of  the  deformatlonal  sequence. 
This  level  of  stress  may  be  termed  the  dynamic  yield  stress  of  the  metal  in 
plane  stress  at  constant  strain  rate.  For  this  case  the  dynamic  strength  of 
the  material  has  a  simple  algebraic  expression  as  a  function  of  strain  rate. 

^  _  A  ^  A _ 

dynamic  yield 

This  expression  is  valid  only  for  strain  rates  such  that 

€  °  >  t 
P  T 

'These  values  of  the  dynamic  yield  strength  for  a  titantium  ed.loy  Ti  2.5 
Sn  5  A1  are  plotted  in  fig,  2.5  along  with  the  plateau  stresses  reached  in 
the  Frank- Read  model  (case  d).  Also  plotted  are  the  experimental  values  for 
the  yield  stress  at  constant  strain  rate'  The  constant  dislocation  model 

gives  dynamic  yield  strengths  which  are  too  high.  The  regeneration  model  gives 
reasonable  agreement  with  the  experimental  data. 
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Section  3 


FRACTURE 

In  the  perforation  of  metallic  plates  hy  projectiles,  it  is  rather  oibvious 
that  the  plate  will  fracture  if  perforation  is  to  occur.  To  determine  when 
a  metal  specimen  will  break  under  rapidly  applied  short-dTiration  loads,  some 
relationship  between  fracture  and  stress  must  be  set  forth. 

It  is  obvious  that  some  combination  of  finite  stresses  will  cause  a  solid 
specimen  to  break  up  since  the  solid  has  only  a  finite  bonding  energy.  It 
is  observed  however  that  the  work  done  to  break  up  a  solid  is  much  less  than 
the  bonding  energy. 

In  1920  A. A.  Griffith'  '  advanced  the  postulate  that  the  strength  of 
a  solid  material  could  be  explained  by  the  presence  and  subsequent  growth  of 
cracks  within  the  specimen.  (That  this  theory  is  fundamentally  correct  has 
been  long  since  verified  experimentally).  Griffith's  postulate  was  that 
the  crack  would  grow  if  the  rate  of  release  of  elastic  strain  energy 
STirroTonding  the  crack  was  greater  than  the  increase  in  the  surface  energy. 

It  is  observed  in  the  fracture  of  metal  that  the  crack  will  grow  under 
applied  stresses  if  the  stress  is  large  enough.  The  stresses  in  the  region 
near  the  crack  are  very  large.  Even  for  so-called  brittle  metals  substantial 
plastic  flow  is  noted  next  to  the  growing  crack. 

(S9) 

Thonqpson  et  al.  have  qualitatively  described  the  fracture  process 
for  polycrystalline  copper.  This  description  might  well  be  an  excellent 
starting  point  for  a  thorough  discussion.  In  their  stressed  specimen,  the 
first  cracks  developed  in  the  glide  bands  that  formed.  The  length  of  these 
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cracks  then  Increased  ^accompanied  by  plastic  deformation  in  the  neighborhood 

of  the  crack.  Then  several  cracks  that  were  formed  Joined  together.  The 

entire  process  of  fracture  under  applied  stresses  may  then  be  broken  Into 

(20) 

3  distinct  phases; 

Phase  I  -  Formation  of  the  crack 
Phase  II  -  Growth  of  the  crack 

Phase  III  -  Joining  of  the  cracks  to  cause  macroscopic  fracture 

Phase  I  or  the  formation  of  the  cracks  themselves  within  a  good  crystal 
matrix  has  not  been  overly  considered  in  the  theory  of  fractures.  In  the 
first  place, the  normal  metallic  specimen  already  has  the  microcracks  in  it. 
Second,  the  formation  of  the  mlcrocracks  is  tied  up  with  plastic  deformation 
on  the  crystalline  level,  a  subject  which  is  not  well  understood. 


Almost  all  subsequent  fracture  theories  of  solid  specimens  have  followed 

( eo 

this  basic  model  and  have  concentrated  on  phase  II.  Sneddon'  '  and  Westergaard^ 

(21) 

calculated  the  elastic  strain  energy  for  various  shaped  cracks,  Irwin' 
post'ulated  that  the  rate  of  energy  given  up  by  the  solid  comes  from  plastic 


strains  near  the  crack. 


Rather  than  try  energy  approaches  to  this  problem  of  the  growing  crack, 
we  shall  ask  one  simple  question.  Given  a  crack,  where  does  it  get  the  empty 
space  necessary  to  grow?  This  question  may  be  answered  in  two  ways.  Flrst^ 
there  is  the  possibility  of  straining  or  displacing  the  walls  outward.  However, 
if  the  stress  is  relieved,  the  volume  of  the  crack  will  return  to  its  original 
position.  To  extend  the  volume  covered  by  the  crack,  either  bonds  must  be 
broken  at  the  edge  of  the  crack  or  dislocations  must  enter  into  the  crack  from 
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the  surroxinding  matrix.  To  hreak  a  bond  requires  a  great  deal  more  energy 
than  to  move  a  dislocation. 

The  model  that  vill  be  developed  formally  in  the  next  section  then  assumes 
phase  I  conq>lete>  that  the  growth  of  a  crack  is  caused  by  dislocations  moving 
into  it  and  that  the  cracks  will"  grow  into  and  overlay  each  other. 

The  dislocations  however  have  a  direct  role  to  the  plastic  strain  of 
the  crack  region  and  to  the  specimen  as  a  whole.  Through  this  feat\ire,  the 
relationship  between  fracture  and  external  forces  is  sought. 

3.1  Mathematical  Model  for  Fracture  under  Transient  Loading 

Consider  an  elementary  volume  of  some  metal.  This  elementary  volume 
contains  initially  mostly  little  volumes  of  perfect  crystalline  material  and 
some  imperfections.  These  imperfections  fall  into  several  categories. 

<•  1.  Interstitial  atoms  and  vacancies 

2,  Intracrystalline  dislocations 

3,  Intercrystalline  dislocations  and  grain  boundaries 

4,  Foreign  or  solute  atoms 

In  addition  to  these  atomic  level  imperfections  there  are  imperfections 
on  a  larger  scale. 

This  elementary  volxime  will  contain  many  individual  crystals  of  the 
metal.  These  crystals  will  have  a  random  orientation  with  respect  to  any 
set  of  Cartesian  axes.  (In  many  cases,  however,  this  may  not  necessarily 
be  so.  The  following  argvments  would  then  have  to  be  altered  to  acco\mt  for 
this  anlsotrojy).  The  main  result  of  this  assunqptlon  is  that  local  variations 
of  Imperfection  movement  will  be  cancelled  out  by  averaging, 
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Let  the  elementary  volume  be  placed  vinder  a  \mlform  stress  from  external 
forces.  The  stress  field  within  this  volume  will  he  con^sed  of  the  external 
stress  and  the  local  stresses  due  to  the  various  Imperfections. 

In  an  elementary  voliane  denote  the  volume  of  one  of  the  macro-holes  hy 

V  •  Two  factors  must  now  he  estahllshed.  First >  given  the  rate  of  growth 
m 

of  a  single  macro- hole,  the  rate  of  growth  of  the  total  volume  V  of  the  macro - 
hole  must  he  computed.  Second,  the  rate  of  growth  of  a  single  macro- hole 
must  he  estahllshed. 

Assume  that  at  some  time  the  volume  of  metal  is  divided  into  two  volximes, 
a  volume  of  remaining  good  crystal  V^-V  and  the  volume  V  of  the  holes.  The 
volume  associated  with  each  elementary  hole  grows  at  some  rate. 

=  k(t)  (3.1) 

The  rate  at  which  the  elementary  hole  grows  is  dependent  on  external  variables 
as  much  as  stress  and  temperature  and  will,  in  general,  vary  with  time. 

The  increment  in  the  volinne  occupied  hy  all  the  holes  is  proportional 
to  the  remaining  vol\ime  of  good  crystal  and  to  the  rate  at  which  an  elementary 
hole  Increases 

dV  =  (V^-V)  dVjj  =  k  dt  (3.2) 

or  expressed  in  terms  of  rates 

dY 

f  ■  (3.3) 

Integrating  the  equation  for  the  increment  in  total  volume  occupied 
hy  holes 
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v= 


1  -  e 


(3.U) 


or 


-I 


(3.5) 


The  next  step  is  to  determine  an  expression  for  V^.  This  is  the  portion 
of  the  problem  that  will  have  many  possible  variants,  i.e.  dependence  on  crystal 
structiire,  stress  i)attem,  temperatvire,  etc.  An  outline  of  the  procedure 
used  in  determining  this  rate  will  be  given  and  only  one  of  the  many  possibilities 
will  be  used  in  further  calculations. 


The  hole  with  volume  has  some  shape.  It  lies  wholly  within  a  region 
of  fairly  good  crystal  containing  several  types  of  atomic-level  imperfections 
(Fig.  3.1).  Under  a  uniform  external  stress  field,  some  of  these  dislocations  will 
move  toward  the  hole.  When  they  intersect  the  hole,  they  will  cause  the  hole 
to  increase  in  size.  (A  dislocation  or  a  vacancy  is  Just  a  very  small  hole 
or  en^jty  space.  What  is  happening  here  is  that  the  volume  of  very  small  holes 
within  the  matrix  of  the  good  crystal  are  being  transferred  to  the  large  hole). 

To  simplify  further  discussion,  the  only  imperfection  to  be  considered  will 
be  edge  dislocations,  which  is  the  dominant  mode  of  inelastic  behavior.  The 
author  has  previously  characterized  the  behavior  of  these  imperfections  under 

Cl3 

transient  external  loads  ,  Because  of  the  initial  assun^ion  of  random 
orientation,  we  need  only  to  consider  the  average  motion  of  the  dislocations. 

The  growth  of  the  micro-holes  requires  that  dislocations  move  through 
the  body  or  matrix.  The  dislocations  will  move  in  patterns  relative  to  the 
local  stress  field.  The  stress  in  the  matrix  is  coiiq>osed  of  stresses  from 
two  sources,  the  applied  stress  and  the  stress  concentration  arooind  the  holes, 
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Dislocation 


Figure  3*1  -  Initial  Configuration  of  Microhole 
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If  the  dislocations  have  to  move  for  any  appreciable  distance  within  the  crystal 
before  striking  a  hole^  the  majority  of  its  motion  will  be  determined  by  the 
applied  field.  We  shall  assume  that  is  the  case. 

The  geometry  of  the  stress  pattern  determines  the  directions  from  which 
the  dislocations  approach  the  holes.  For  example,  in  uniform  compression, 
the  dislocations  will  approach  the  hole  in  a  45“  cone  above  or  around  axis 
of  largest  principal  stress'^'CPig,  3*2).  For  simple  shear,  edge  dislocations 
will  come  only  from  one  direction  (parallel  to  the  shear  in  traction)  Fig.  3»3* 


The  hole  will  grow  at  a  rate  proportional  to  the  surface  area  exposed 
to  oncoming  dislocation.  The  hole  will  grow  at  a  rate  proportional  to  the 
number  of  dislocations  that  cross  this  area  per  unit  time.  It  is  also  pro¬ 
portional  to  the  volume  increment  of  each  single  dislocation 


(3.6) 


Ag  =  surface  area  of  hole  exposed  to  oncoming  dislocations 
N  =  nxanber  of  dislocation  per  unit  area 


AV  is  voltune  associated  with  a  single  dislocation 


If  a  dislocation  has  a  velocity  Vjj  it  will  reach  the  hole  in  unit  time 
if  it  is  close  enough.  This  requirement  defines  a  length  X  such  that  all  of 
the  dislocations  in  the  volume  of  good  material 

(3.7) 

will  reach  the  macro-hole. 


Thus  the  rate  of  dislocations  crossliog  into  the  hole  may  be  written  as 


dt  ~T~ 


(3.8) 
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New  Line 


Original  Hole 


/  V^. 


Dislocation 


Boundary  of  the  cone  in 
which  lie  the  dislocations 
that  will  intersect  the  hole 


Figure  3.2 

Growth  of  the  Microhole  \inder  TriaxisQ.  Compressive  Stress 
(Spherical  Growth) 
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Figure  3.3 


Growth  of  the  Microhole  under  Single  Shear 
(Rod-like  Growth) 
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For  the  compression  case,  the  expansion  of  Is  spherical 


^  _  llTrlr  T*^  A 
- r  A  -5;- 


(3.9) 


dr  V 


(3.10) 


t 

^  -  *'0  =  /  ^2  -r- 

o 


>^2  -  ‘‘1'^  ^dls 


The  volume  of  the  hole  at  any  later  time  Is 


/i4f  t  V. 

\  =  T  ^  ^2 


VdN 

-T- 


with  V  the  initial  volvime  of  the  hole, 
tn 


(3,11) 


For  simple  shear  case  the  growth  Is  rod-like.  The  linear  rate  of  growth  is 


dx  , 

cT^  “  “T" 


The  volume  of  the  hole  at  any  later  time  is 


(3.12) 


V  = 
m  s 


dt  +  V_ 


(3.13) 


Inserting  these  expressions  into  the  equation  for  total  volume  of  holes 


r  ^  ^ 

-vrf.at)3 

V  =  V  1  -  e 
o  L 


spherical  or  cor^jressive 


V  =  LI 


Rod-like  or  simple  shear 
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Xf  the  stress  is  veiy  hlgh>  Vjj  Is  very  nearly  constant  and  the  integrals 
may  be  evaluated  explicity 


-V 

V  =  V^(l  -  e  ^  ) 


Conipressive  stress 


with  ki.  = 


V  =  V  (1 
o 

V  °N 

^  ''d  " 


-k  t 
e  5  ) 


Simple  Shear 


’‘5  ■  '‘3  X 


A  schematic  picture  of  the  rate  of  growth  of  the  totsuL  volume  occupied 
by  the  micro-holes  is  shown  in  fig.  3.^. 

These  two  curves  may  cross  each  other.  It  is  quite  likely  that  the  two 
coefficients  k|^  and  k^  are  different  functions  of  temperature.  ThuSj  at  one 
temperature  of  one  of  the  mechanisms  may  grow  much  faster  than  the  other# 
while  at  a  different  temperature  the  order  may  he  reversed. 

3.2  Relation  of  the  Fracture  Criteria  to  Plastic  Strain 

The  postulate  for  fracture  of  the  elementary  volxune  is  very  simple. 

If  the  volume  occupied  by  the  assemblage  of  macro-holes  is  greater  than  a 
certain  percentage  of  the  total  volume  the  structures  can  no  longer  hold  to¬ 
gether  and  the  volume  falls  apart  —  or  in  other  words#  fracttires.  This 
critical  volume  ratio  may  be  different  for  different  crystal  structures. 

This  critical  volume  concept  holds  only  In  the  region  of  the  specimen 
wherein  the  shear  stresses  are  high.  The  growth  of  the  elementary  volumes 
is  dependent  on  the' level  of  the  resolved  shear  stress.  In  any  large  specimen 
these  shear  stresses  are  not  uniform  throughout.  At  the  beginning  of  the 
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Figure 

Stress  vs,  TJ rae-tc-Fi’acture  at  Constant  Stress 
Coii5)ared  to  Plastic  Strain-Rate  at  Constant  Stress 
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deformation  process  when  all  the  stresses  are  elastic^  this  stress  pattern 
Is  not  iinlform.  As  plastic  deformation  of  the  specimen  proceeds  this  non- 
unlformlty  of  stress  distribution  and  concentration  can  become  altered.  In 
particular,  the  high  shears  can  localize  In  or  near  the  area  of  greatest  plastic 
flow  (for  a  numerical  example  of  this,  see  ref.  12).  The  shear  stress  pattern 
must  be  solved  as  a  part  of  the  boundary  value  problem  involving  the  elastic 
strain  and  flow  pattern. 

The  critical  volume  for  fracture  may  be  quite  different  for  different 
types  of  crack  growth,  for  example,  the  crack  that  grows  in  a  disk  shape  will 
require  only  a  small  volinne  before  fracture  occurs. 

This  critical  volume  will  be  epsentially  the  volume  of  holes  in  a  region 
wherein  the  shear  stress  is  uniform  and  will  not  pertain  to  the  volume  of  holes 
in  the  entire  specimen. 

Denote  by  V  the  ratio  of  V  to  V  when  fracture  occurs.  At  fracture 
“  c  o 

-V 

Tr  n  m 

V  =  1  -  e 
c 


=  1  -  e 


.  dV 

7  ^  dt 

•*0  dt 


dV 


m 


Inserting  the  expression  for 


V  = 
c 


Vd 


n 

Ndt) 


1  -  e 


(3.14) 


where 


n  =  1  for  simple  shear 


n  =  3  for  biaxial  con^ression 
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Solve  this  eqtaation  for  the  Integral 
t  -l/n 

=  [i°(l-Vc)]  (3.15) 

This  Integral  expression  may  he  compared  to  the  expression  for  toteil 
plastic  strain  of  the  previous  chapter, 
t 

e  =  h  /  dt  (3.16) 

^  o 

where  h  is  the  Burger's  vector  defined  in  Sec.  2. 

Thus 

"1/^ 

=  Xir  tin(l-v^)  ]  (3,17) 

fracture  ^ 

Thus  at  fracture  a  certain  amount  of  plastic  strain  has  occured  which 
is  an  indicator  of  fracture.  The  magnitude  of  this  critical  strain  is  one  measure 
of  whether  the  metal  is  ductile  or  brittle. 

Since  the  plastic  strain  rate  for  simple  compression  or  tensile  stress 
specimens  is  a  function  of  the  applied  stress,  experiments  carried  out  at 
constant  stress  sho\ild  give  some  indication  as  to  the  validity  of  the  theory. 
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Stress*^  ksl 


At  some  constant  stress  level.*  the  plastic  strain  rate  Is  constant  And 


using  the  form  for  this  strain  rate  as  In  Sactlon  2. 


=  V_,Nbe 

P  D 


-A/a 


(3.18) 


where  a  is  the  applied  stress. 


Integrate  the  plastic  strain  from  zero  to  the  critical  strain 


e  °  =  e  (a)  t„ 
p  p'  '  f 


(3.19) 


where  t^  is  the  time  at  which  fracture  occurs. 
Thus 


tj.  = 


•(3;20) 


A/a 

e 


(3.21) 


This  time-to- fracture  curve  is  very  similar  to  the  data  on  aluminum  and 

(32)05) 

the  noble  metals  measured  by  Zhurkov  .  The  results  of  the  theory  and 

the  experiment  are  shown  in  Fig.  3.5. 

The  tlme-to-fracture  of  metallic  specimen  under  a  given  stress  leading 

can  be  related  to  the  plastic  strain  in  the  region  of  the  fracture  under  the 

same  stress  history.  This  criteria  will  be  used  to  determine  the  penetration. 

/f  o  c 

This  criterion  is  /  ~  G 

o  P  P 

To  use  this  crltei^ion  in  a  boundary  value  problem#  thfe  deformatlonal  process 

must  be  solved  for  the  dynamic  elastic  plastic  case  and  the  fracture  will  occiir 

at  the  point  or  surface  in  the  material  whore  first  the  plastic  strain  exceeds 
Its  maximum  vaXii^« 
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Section  k 


TRAMSIEWT  WAVE  ANALYSIS 

In  this  section  ve  shall  Investigate  the  propagation  of  stress  waves 
through  an  elastic  medium.  First  we  consider  an  axled-ly  symmetric  half- space 
\inder  an  axhltrary  axially- symmetric  load.  The  general  Integral  solution  for 
this  problem  is  then  examined  for  the  specific  case  of  a  \mlt  load  over  a 
unit  circle,  and  expressions  are  derived  for  the  stresses  along  the  axis  under 
this  load,  as  functions  of  time  and  depth.  A  numerical  integration  is  programed 
and  carried  out  on  the  IBM  l620  to  find  the  values  of  these  stresses. 


A  second  program  is  then  developed  which,  from  the  numerical  values 
Just  found,  computes  the  stresses  as  the  waves  are  reflected  through  a  plate. 
This  con^nitatlon  is  performed  for  two  cases,  a  Pe  plate  and  a  T1  plate.  The 
results  are  presented  as  contour  drawings,  which  show  the  time,  depth,  and 
duration  of  maximum  stress  cor-jonents  in  the  radleil  and  axial  directions;  the 
tangential  stress  along  the  axis  will  have  been  shown  to  be  identically  equal 
to  the  radial,  and  the  shear  to  be  identically  zero. 


4.1  Equations  of  Motion 

Stress-wave  propagation  is  defined  by  the  equations  of  motion 


rr  ,  rz  ^  rr  yy  _  ,, 

+  -  +  - - — ■  -  p 


Sr 

St 


Sz 


3t‘ 


T  Sa 

£2  +  =  ij. 


Sr 


Sz 


S^u 


St' 


(4.1) 


(4.2) 
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wherfe  ^zz*^66  ^rz  radial,  axial  (v«rtlcaJ.)*  teuagentlal, 

and  shear  stzHfeiBSes  glvdn  by  the  stress-strain  relations 


hu 

u 

a 

rr 

<X^  2G) 

_ r 

Sr 

♦  x(— 

r 

Su.. 

Su_ 

^ee 

<X-f2G5 

-i:  ♦ 

r 

Sr  Sz 

u 

^zz  * 

(X'f2G) 

— ^  + 
Sz 

r 

<U.5) 

G( - -  4 

^z 

Sr 

with  and  denoting  the  radial  and  axial  components  of  displacement,  and 
/i,  G  the  density,  uniaxial  strain  modulus,  and  shear  modulus  of  the  medium. 

h.2  Boundary  Conditions 

The  condition  of  axial  symmetry  is  assumed  in  the  hbove  relations;  the 
further  boundary  conditions  for  this  problem  are  that  the  surface  ii-  =  0  (taking 
the  axis  to  be  vertical)  is  free  from  shear  and  that  the  pressure  history  on 
the  STiTfac©  is  giv^n  by  some  fvinction  f.  Then,  considering  the  stresses  as 
functions  of  r,  z,  and  t,  we  have  for  our  boundary  conditions: 
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The  variables  in  the  problem  can  be  reduced  to  dimensionless  fonn  by 

introducing  a  characteristic  pressure  p  and  characteristic  distance  R.  The 

characteristic  distance  will  be  taken  as  the  radius  of  the  biillet.  The 

f  iW 

characteristic  pressure  is  the  normal  pressure  generated  at  impact'  ' 


where 


p  ^  MiCxMaCg  Y 
^  /LtiCi  +  /iaca  S 


M1M2  are  densities  of  the  plate  and  bullet 

C1C2  are  longitudinal  sonic  velocities  of  the  plate  and  bullet 
Vg  is  the  in^jact  velocity 


In  terms  of  the  dimensionless  variables  defined  above,  the  equations 
of  motion  become: 


3\iatlons  in  Dimensionless  Form 
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♦  i  — £  -  ^  4.  - i  +  — a  =  /s^  — ^ 

P  ^P  p^  ^C^P 

■U-  ,  ^u,.  N  ^,u_  u„  c 


(4.9) 


d^u-  1  du»  N  du^  u^  d^u.  d‘u„ 

^  ^  ^  ^  (4.10) 

dp  p  op  OQ  op  p  o^  or 

And  thfe  "boundary  conditions  nte:- 
Of^l^ipy^tr)  =  -ipip^r) 

Tp^(PrO,T)  =  0 

4*5  Qenoral  Integral  Solution 

The  Hankel -Laplace  transform  solutions  for  arbitrary  f  are  given 

-  4i  I  /"  ^  -  4^)e-'"«]  J.(P?)  de 

Br  o  L  J  (4.1L) 

+/"-P  [(m=  .  -e'mne-"«]  Jo(p?)<i?  ds 

°CC  “  ■  ■ST  ^  J"”?  *  ^~'>  ®’'”^  mne""^]  Jo(p5)  d^ds  (4i2.) 


mne"^^  Jo(p|)  d^ds  (1^X20 


PP^’^00  "2^ 


k  L  / 


{e=-o=-2)s=)(e=  +  ■e5s!)e-”^-5 


jo(pe)deds 


(413) 


PC  sr- g 


(2C^  +  (n^  +  Ji(pC)  d^ds 

-  •  /V. 


(4.a4y) 


where 


2 

g  =  (4^  +  -^-1-)  -4^mn 


?  #  /"  /"  e"®''’  0(p,t)p  Jo(p4)  dpdT 


MRD  DIVISION 

OENERAI*  AMERICAN  TRANSPORTATION  CORPORATION 


h,6  Specific  Infgral  Solution 


If  the  presBtire  function  Is  a  imlt  pitessure  (p)  applied  Instantaneously 
and  held  constant  over  a  circular  area  of  unit  radius  (r),  we  have 

to  ' 

?  =  /  6“®”^  /  P  Jo<p4)  dpdT 

o  o 

o 

^  a  £A(,i) 

CS 


We  seek  the  stresses  generated  aloxig  the  axls^rp=0)  hy  this 'pressure. 
Since  J^(o)  =  0,  we  have  from 


s  0 

In  evaluating  cr  we  note  the  factor  *^^(^4). 

H  H  H 

Since  lim  J-(p^)  =  0,  we  apply  L’ Hospital's  rule; 

p+0 

p-»o  P  p-»o  ^  P 


in  the  first  term  of 


=  llm  e  JiCpO 

p^o  ^ 


«  aim  e 

P*o 

=  Jo(p^) 


wherein  we  retain  this  Bessel  function  for  later  use. 
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Thus  we  have 


(O.C.T)  =  f  !£  I-  -  |-(e=  +  -e  nme-"? 


i.  ^ 

t  -  £^)<4=  +  Jo(pO  dids 


4r  J'”l|  (?^-0^-2)s^)(e  +  £^)e-”'?-e=  mne-"«  Jo(p?)  d^ds 


And  hence  at  p  =  0 


^pp  ”  ^<j)4>  ^^^pp  ^  *^00^ 


Then  recalling  that 


if  = 


we  must  evaluate  the  two  integrals 


2^pp-2a^^=  ^  f^  — 


-Jo(pOJiU)deds 


^  mne- 


/  I-/ 


Jo<pe)  Ji(Od|ds  (Mr) 


We  can  eval\aate  the  three  integrals 


1  /  £I 

^  Br  ®  =  * 


Jo{pO  JiCOd^ds 


1  ^ST  »  (^2_f^2_2)s2)C^^  +  .  I 

l2  =  /  s —  /  - - - g Jo(pO  Ji(Od|ds/ 

Br  o  ° 

,  ST  W  *2  f. 

I±  =  J  I—  /  Jo(pe)  Ji(e)deds 


o  e 
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and  \rm  have 


('O^  C>'^)  =  la  “  Ii 


CJpp(0,C,T)  ^  =  ^l2  -  I3] 


(4.19) 


4.7  Analytic  Reduction  of  Integral  Solution 

We  shall  deal  simultaneously  with  the  three  Integrals  of  (4,l8).  Let 
us  define  a  new  variable  x  as  follows:  on  the  Bromwich  contour  s  =  e  +  16. 
For  6  <0  let  ^=sSXj  for  6^  let  ^  =^-sx.  This  gives  d4=sdx  and  d^=-sdx 
respectively.  Since  ^_^Re(x)  =  0f  we  can  consider  the  path  of  integration 

for  X  to  be  the  imaginary  axis,  and  our  choice  for  the  sign  makes  the  limits 
0  and  loo  for  both  integrals.  Then  taking  Ii  for  example. 


ds  . . 

a-TTl  e-l«  0  (x^  +  jSVs)^  -  v^tx^'+l)  (x^+^ 


[e"®^  Jo(psx)  Ji(sx)  dx]  +  ds} 

r\ _ j 


(rl~  [ _ (x^  +  /3V2)^  1  [-^-sC  A^+1  j 

®  (x^  +  /3^/2)  ^  -  x^  ■/(T‘^+1)  (x^+jS^) 


[Jo(-psx)  Ji(-sx) (-dx) ]} 


SlhC)»  Jo(psx)  and  Ji(sx)  are  even  and  odd  functions  respectively,  the  minus 
signs  will  drop  out.  Reasoning  similarly  for  Ig  and  142,3 


BT.„  „-sC  V^c^+oT* 


■*"  27rl 


P(x) 


k  Jo(psx)  Ji(sx)  dx  (4.20) 
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where 


fi(x)  =  (x^  ♦  I”) 


f2(x)  = 

1 1 

r*) 

f3(x)  = 

ri 

F  (x)  = 

tti  =  02 

F  (x)  =  (x®  +  ■^)  -x^\/(x^+l)  (x^+jB^) 


ttg  =  i3 


Now  we  make  the  substitution  Jj^(sx)  =  ^(H^^^(sx)+H^^^(sx))  and  let  x=lv 
dx=idv ; 

T  1  r  -ST,,  “  S^(v)  -st/a^=-v=('H<‘‘)(lsv)+Ht^)(lsv)] 

^  'SttT  J  ®  J  TTHTT-  ®  1-^^ - 5 — * - 


Jo (ipsv) dv 


where 


gi(v)  =  (| - 

g2(v)  =  (2-/3^-v^)(| - v^) 

g3(v)  =  -v^-/(l-v^)  (j3^-v^) 

G(v)  =  (|^  -  v=.)^  +  vV(l-v2)(j3=-v2) 
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Splitting  the  Integral  Into  tvo  terms  and  changing  the  second  to  the 
negative  half  of  the  axle  gives 


ds  [i  /  e  ^  ^  ~  Jo{lpsv)  dv 


■2  J  gTv] 

o  '  ‘ 


2  J 

o  ' 


(isv)  Jo(lpsv)  dv 


Noting  that  H^^^(lx)  =  -  ^  k^(-x),  and  that  I^(x)=J^(lx), 


we  have 


/  e®'"  ds  i  / 


«  gi^(v)  -siya^^-v‘ 


Ki(sv)  lo(psv)  dv 


-00  g,  (v)  -sCVo,  ^-V^ 

^  G^v)""  ®  ^  Ki(-sv)  lo(psv)  dv 


which  by  Bateman's  theorem  on  conjugate  integrals'  '  reduced  to 

,  CO  g,  (v)  -sCva,  ^-v^ 

-  ski  IT  ^  GT^  ®  Ki(sv)  lo(psv)  dv 

Br  o''  4  ioo 


Complex 

V-plane 
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Considering  the  Inner  integration  as  integration  in  the  convex  v  plane> 
we  note  that  there  are  no  singularities  in  the  right  half;  let  us  consider 
the  integration  over  the  path  indicated  above,  from  0  to  oo ,  thence  to  ioo 
and  hack  to  zero  along  the  right  side  of  the  imaginary  axis.  The  Integral 
over  the  closed  path  is  zero,  and  since  the  Integrand  vanishes  on  the  portion 

O  A 

of  the  path  from  oo  to  ioj ,  we  have  /*+  /.  =  0  or  /**=  /  **  .  Hence 

o  !•“  o  o 

can  he  written  as 


_  If  ST  ,  21  r 

=  ■  -^^rT  ®  ^  ^  o 

or  inverting  the  order  of  integration 


1«  Siri^)  -sC/a^^-v^ 


Ki(sv)  lo(psv)  dv 


^  ai  1 

—  J  TJtVr  ■SFT 

o  '  ' 


sCT-C-Za.  ^-v^) 

/  e  Ki(sv)  lo(psv)  ds 

Br 


We  now  let  v  =  lx,  s  =  -iy,  yielding 

00  f,  (x)  oo  -lyT), 

dx  /  e  ^  Ki(xy)  Io<pxy)  dy 


where  7]^=t-  4  and  and  F(x)  are  defined  as  in  (4£0')« 

But  =  cos  yn-i  sin  y?l^o  the  real  part  of  the  integral  is 

-  00  f|^(x)  p  « 

Re(lk)  ==  w  f(x')  ^  tt  ^  Ki(xy)  lo(pxy)  dy 


Now  hy  a  theorem  of  Sneddon 


(10) 


^  /  sin  yr)  Ki(xy)  lo(pxy)  dy  =  /  Ji(xy)  Jo(pxy)  dy 
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2 

forT)  »‘0.  This  sets  an  upper  Halt  of  x  for  the  outer  Integral^ 


fiNd  p 


^  fj^<x)  •  -s<t-Cv4j^^+x2) 

Re(lj^)  «  ^  /  TOT"  ^  ®  Ji(sx)  ds 

The  inner  integral  is  a  standard  form;  evaluating  it  gives  us  the  desired 
solution 


Re(li,)  = 


TT  •*  u  P(x) 

O 


)  ^  +  X^  -  (T-CVttj^^+X^) 

x'V(T-Cv^ar^x^^^' + 


dx  (4.B(fe) 


U.8  Kumerlcsil  Integration 

The  integrals  were  evaluated  numerically  on  the  IBM  l620.  Values 

of  T:  ran  from  .1  to  3»0  at  intervals  of  .1,  and  values  of  C  were  from  .1  to 
T-*lf  also  at  .1  Intervals.  The  integration  was  also  performed  for  C  =  'T”-05j 
T  .02,>t-.01  and  T-.005  >  to  give  an  approximation  to  the  shape  of  the 
peak  and  the  value  at  the  wave  front.  The  parameter was  taken  to  he  3- 

2 

For  each  point  in  the  T-^mesh#  the  upper  limits-**^  -1  and  were 

computed  and  designated  as  ULl  and  UL2  respectively  (UL2  was  set  eqvial  to 

zero  for  The  program  then  defined  DXL  =  and  DX2  =  The 

values  of  the  Integrand  were  found  at  the  points  x  =  n»EQQ.  (n=l,  .,.,20)  for 
the  first  two  Integrals  and  x  =  n»nX2  for  the  third.  The  Integrals  were 
then  approximated  hy  Sls^son's  rule>  and  the  stresses  found  by  adding  the 
appropriate  Integrals. 
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The  results  of  this  con^nitatlon  are  shown  In  Figure  4.1,  for  selected 
values  of  T,  with  o^^and  plotted  as  functions  of 

4,9  Reflection  of  Stress  Waves  ■’ 

The  radial  and  aocial  stress  components  have  now  been  computed  for  waves 
propagating  in  a  half- space.  To  apply  these  result's  to  a  plate  of  finite 
thickness,  a  computer  program  was  developed  to  find  the  stresses  reflected 
in  the  plate  from  the  stresses  mentioned  above.  That  is,  the  half- space 
stress  previously  computed  were  used  as  inputs  for  this  program  which  for 
each  point  in  the  time-depth  mesh  calcoilates  the  sum  of  the  stresses  from 
all  reflected  stress  xreives  which  have  passed  the  point,  after  multiplying 
each  of  these  stresses  by  the  product  of  the  appropriate  reflection  coefficients 
to  accoxmt  for  the  reflections  undergone  by  the  wave  in  question.  Since  the 
stresses  are  calculated  initially  in  terms  of  a  unit  pressure  input,  the 
final  step  in  the  computation  is  to  multiply  the  stress  sum  by  the  coefficient 
of  transmission  from  the  projectile  into  the  plate.  The  materials  considered 
here  are  Fe  for  the  projectile  and  Fe  and  Ti  for  the  plate.  The  transmission 
and  reflection  coefficients  for  these  media  are  as  follows^  ^)[22) 


o 


Fe 


Ti 


V 
+> 

+» 

V^Siicuum 
a 


i 


Incident  Wave  in 


Ti - 

- 

T  =  1 

T  =  .735 

R  =  0 

R  =  ,265 

T  =  .735 

T  =  1 

R  =-.265 

R  =  0 

T  =  0 

T  =  0 

R  =  -1 

R  =  -1 
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where  transmission  into  a  vacrium  corresponds  to  our  assusqptlon  that  the  back 
of  the  plate  is  free.  The  dimensionless  id.ate  thickness  C  is  taken  as  0.8 
for  each  plate/  and  stresses  are  cooqputed  for  O  2  at  intervals  of  .10. 

In  the  case  of  the  Fe  plate  (Figure  4.2),  only  the  incident  wave  from  the 


projectile  and  the  first  reflection  from  the  back  of  the  plate  appear,  since 


c 

Figure  4.2  -  Reflection  of  Stress  Waves  in  Fe  Plate 


the  latter,  upon  reaching  the  front  of  the  plate,  is  entirely  transmitted 
into  the  projectile.  Let  us  consider  for  example  the  social  stress  at  point 
A.  Since  this  point  lies  behind  both  the  Incident  wave  front  and  the  reflected 
wave  front,  there  will  be  two  terms;  first,  the  direct  stress  at  A;  second, 
the  product  of  the  Fe-^Vac  reflection  coefficient  by  the  stress  at  the 
reflection  A*  of  A  about  the  rear  of  the  plate.  The  sum  is  then  multiplied 
by  the  projectile  -*•  plate  transmission  coefficient: 

a^^(2.0,0.2)  =  [a'^^(2.0,0.2)  +  (-l) *  a« ^^(2.0,1.4) ]  *  1 
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vhere  a denotes  the  stress  coopted  am  the  half- space* 

For  the  corresponding  point  T1  plate,  there  are  three  waves  to  consider 
(Flgiire  4,3),  since  the  reflection  from  the  hack  of  the  plate  is  not  entirely 
transmitted  into  the  projectile  when  it  reaches  the  front  surface,  hut  is 
partly  reflected  hack  into  the  plate. 


Thus  we  have  in  this  case 

^^{2.0,0.2)=[a'^^(2,0,0.2)  +  (-l).  a«  ^^(2.q,  1 . 4)  +  (-lMR^^^pg)  .  o' ^^X2. 0,1 .8) 


where  T  and  R  are  the  transmission  and  reflection  coefficients  as  indicated 
hy  the  subscripts. 


The  resxalts  of  this  computation  are  shown  as  stress  contours  in  Figures 
4.4  and  4.5,  plotted  in  theT-?  plane.  The  heavy  diagonal  lines  in  these 
two  figures  show  the  fronts  of  the  incident  and  reflected  waves.  The  shaded 
areas  indicate  regions  of  tensile  stress. 
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4.10  Wave  Analysis  for  Laminated  Plate 


The  anaJ.y8ls  discussed  above  Is  extended  here  to  provide  a  solution 
for  the  stresses  In  a  laminated  plate;  we  consider  plates  consisting  of  a 
layer  of  Fe  and  a  layiir  of  Tl,  In  both  orders  and  various  relative  thicknesses. 

The  necessity  for  a  somewhat  simplified  approach  can  be  seen  by  referring 
to  Appendix  A,  where  the  formal  Integral  solution  for  a  laminated  plate  Is 
derived  \inder  a  set  of  rather  artificial  boundary  conditions  chosen  largely 
for  the  relative  sln^jllclty  they  provide  while  still  retaining  some  resemblance 
to  the  physical  slt\iatlon.  Even  under  these  conditions.  It  can  be  seen  that 
the  formal  solution  Is  of  an  Inordinate  degree  of  complexity:  evaluation 
requires  that  8  equations  similar  to  (A. 38)  be  solved,  and  each  of  these, In 
order  to  be  solved,  must  be  further  complicated  by  substitution  of  the  defined 
values  of  the  variables  m^,  n^^,  p^^,  q^,  and  by  the  Insertion  of  an  expression 
for  the  transform  of  the  pressure  history.  Clearly,  the  practical  value  of 
this  formal  solution  Is  somewhat  limited, 

4.11  Method  of  Solution 

As  stated  above,  the  approach  used  Is  the  same  as  that  previously  discussed 
for  a  single  plate.  The  half- space  stresses  computed  there  were  again  used 
as  Input  data  for  determination  of  the  plate  stresses.  The  computation  was 
carried  out  for  six  cases:  the  two  plate  materials  considered  were  again 
Fe  and  Tl,  and  In  each  order  of  the  two,  the  stresses  were  computed  for 
conqxsnent-plate  thickness  ratios  of  1:3#  Itl,  and  3;1* 


■7 
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Numerical  Evaluation 


Detailed  discussion  of  the  program  for  the  second  phase  of  the  problem, 
the  conrposltlon  of  the  stresses,  has  been  reserved  for  this  section  because 
the  same  program  was  used  here  as  for  the  single  plates;  the  program  was 
act\ially  written  for  the  laminated  plates  and  the  single  plates  treated  as 
a  special  case. 


if.  12,1  Input  Data 


There  are  two  types  of  data  for  the  program;  first,  the  half¬ 
space  stresses  computed  In  the  first  phase  of  this  problem;  second,  the 
characteristics  of  the  plate.  They  are  handled  in  this  order;  first,  the 
stress  data  are  read  (in  the  part  of  the  program!  enclosed  on  the  block  diagram 
in  dashed  lines),  consisting  simply  of  a  list  of  either  all  radial  or  all 
axial  stresses.  Then  the  plate  data  are  read.  We  assume  that  the  transmission 

and  reflection  of  the  stress  waves  is  described  by  the  transmission  and 

(22) 

reflection  coefficients  for  plane  waves  as  derived  by  Kolsky' 
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The  R  and  T  read  In  are  simply  the  reflection  and  transmission  coefficients 
Indicated  hy  the  subscripts  0  for  the  projectile,  1  for  the  upper  plate, and 
2  for  the  lower;  e.g.,  Rg^  Is  the  reflection  coefficient  for  a  wave  In  the 
lower  plate  reflecting  frcm  the  upper.  It  Is  noted  that  the  reflection 
coefficient  at  the  bottom  surface  of  the  plate  Is  -1  for  either  plate  materleJ.. 
The  code  number  KHl  Is  3>  5i  It  or  9  for  an  upper  plate  thickness  of  .2,  .4, 

.6,  or  ,8  respectively;  the  total  plate  thickness  Is  .8  In  all  cases. 

U.12.2  feutl^he  of  .Cofl^t'er  Progiton  • 

The  program  for  this  phase  Is  diagramed  In  Figures  4.6  -  4.8. 

After  the  data  are  read  as  described  above,  the  numbers  HI,  H2  and  KH  are 
computed.  The  first  two  of  these  are  the  thicknesses  of  the  top  and  bottom 
plates,  and  are  computed  only  to  be  printed  so  that  the  output  can  be  Identified. 
The  actual  values  of  time,  dejxth,  and  plate  thickness  are  not  used  In  the 
computations  nor  In  the  logical  decisions,  which  can  be  performed  more 
conveniently  by  working  with  the  fixed-point  Indices  corresponding  to  these, 
defined  by 

I  =  lOT  +  1 

J  =  lOZ  +  1 
KHL  =  lOHL  +  1 

The  additional  constant  KH  Is  defined  because  of  the  convenience  of  having 
a  set  of  values  1,  2,  3#  4  corresponding  to  the  vEJ.ues  of  EL, 

The  first  step  In  the  program  after  the  data  and  required  constants 
are  eill  available  Is  to  Initialize  1  =  2.  T  Is  then  computed  (initially  .l) 
and  printed.  Determination  of  UK  corresponds  to  the  statement  that,  for 
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The  R  and  T  z«ad  In  are  8lii:5)ly  the  reflection  and  transmission  coefficients 
Indicated  hy  the  subscripts  0  for  the  projectile,  1  for  the  upper  plate,and 
2  for  the  lower;  e.g.,  is  the  reflection  coefficient  for  a  wave  In  the 
lower  plate  reflecting  from  the  upper.  It  Is  noted  that  the  reflection 
coefficient  at  the  bottom  surface  of  the  plate  Is  -1  for  either  plate  material. 
The  code  number  KHl  is  3>  5$  7>  or  9  for  an  upper  plate  thickness  of  .2,  .1)-, 

.6,  or  ,8  respectively;  the  total,  plate  thickness  is  .8  in  all  cases. 

4.12.2  i&utl'ihe  ctf  .Cbfi^t'er  Program  . *  •  ■ 

The  program  for  this  phase  is  diagramed  in  Figures  4.6  -  4.8. 

After  the  data  are  read  as  described  above,  the  numbers  HI,  H2  and  KH  are 
conqjuted.  The  first  two  of  these  are  the  thicknesses  of  the  top  and  bottom 
plates,  and  are  confuted  only  to  be  printed  so  that  the  output  can  be  identified. 
The  actual  values  of  time,  depth,  and  plate  thickness  are  not  used  in  the 
computations  nor  in  the  logical  decisions,  which  can  be  performed  more 
conveniently  by  working  with  the  fixed-point  indices  corresponding  to  these, 
defined  by 

I  =  lOT  +  1 

J  =  lOZ  +  1 

KHl  =  lOKL  +  1 

The  additional  constant  KH  is  defined  because  of  the  convenience  of  having 
a  set  of  values  1,  2,  3^  ^  corresponding  to  the  values  of  EL, 

The  first  step  in  the  program  after  the  data  and  required  constants 
are  all  available  is  to  initialize  1=2.  T  is  then  computed  (initially  .l) 
and  printed.  Determination  of  NK  corresponds  to  the  statement  that,  for 
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any  T,  Z  will  not  exceed  min  (T,  ,8).  Then  J  Is  set  to  1  (Z=0),  and  J  Is 
con^jared  to  KHl  and  the  d\amny  variable  IDIOT  given  the  value  1  or  2  depending 
on  whether  the  point  being  considered  is  in  the  upper  or  lower  plate  (at 
the  Interface  J  =  KHl(Z=Hl)j  the  stress  is  conqputed  first  for  a  ixjint  in  the 
upper  plate  approaching  the  interface).  The  actvial  stress  computation  for 
points  in  the  upper  and  lower  plates  is  carried  out  in  the  routines  labelled 
and  respectively.  The  sum  of  the  stresses,  S,  is  then  multiplied  by 
TOl,  since  the  half-space  stresses  are  for  a  unit  stress  input,  and  Z  and  S 
are  printed.  The  three  tests  shown  Just  below  this  in  the  lower  right 
comer  of  Figure  4.6  are  to  determine  whether  the  stress  Just  computed  was 
for  the  upper  plate  at  the  interface.  If  so,  the  next  computation  is  for 
the  lower  plate  at  the  same  point.  In  any  other  case,  J  is  tested  to  see  if 
Z  has  reached  its  maximum  value.  If  not,  J  is  incremented  and  the  computation 
is  repeated;  if  so,  I  is  incremented  and  J  reset  to  1.  When  I  is  found  to 
have  reached  its  maximum,  the  computation  has  been  completed  for  the  entire 
plate,  and  the  data  for  the  next  plate  are  read  and  the  program  (except 
for  reading  the  stress  data)  repeated.  It  should  be  noted  that  the  program 
does  not  have  a  logical  end.  After  all  the  plate  data  desired  have  been 
read  and  processed,  the  computer  stops  and  is  then  Instructed  manually  to 
return  to  the  beginning  and  read  the  other  set  of  stresses,  then  go  on  to 
reread  each  set  of  plate  data  and  conqpute  the  other  set  of  plate  stresses, 

4,12.3  Computational  Portion  of  Progiam 

As  has  been  mentioned,  the  composition  of  stresses  is  performed 
in  routines  and  of  the  program,  shown  in  Figures  4.7  and  4.8.  Here 
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Boter 


S=T12 
Index  =  2 
RbT12. 


T 


S*-Si^2*ia.0*Tl2.  gl(l,  J-i-U-KH) 


T;  »;.KH-J+18 


Tt>4.KH-J+18 


Sf6-Rl2 -  R10*T12 • Sl( I, U • KH- J+l8 ) 


K-R-Ii2l 

S«-S+R .  Sl(  I,  Inaex*  (  9-KIIl  )+j) 
Iridex-»Znde::+2 


rcj+i6.Kn 


T;  cT  +  l6 


I  <  J  +  l6 


T  I  i  J  +  l6 

[S^S-R12.T.12  •  K1.0-R21  Sl(l,  J+l6) 


Return 


Figure  4.8  -  DetaiJ  of  Block  Dlsgran 
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we  shall  hrlefly  explain  ,  though  slightly  more  conqjlex.  Is  analogous 


thereto. 


Initially^  the  stress  (radial  or  vertices.)  Is  set  eqvial  to  the  corres¬ 
ponding  half-space  stress;  this  accoiints  for  the  InltleJ.  stress  wave.  Each 
test  In  this  part  of  the  program,  with  the  exception  of  the  test  on  KH,  Is  to 
determine  whether  a  relation  between  1  and  J  (and  hence  between  T  and  Z)  Is 
such  that  the  point  (l,J)  lies  behind  the  front  of  a  iiartlcular  wave.  If  so, 
a  stress  component,  consisting  of  the  product  of  the  half- space  stress  at 
an  appropriate  point  by  the  appropriate  coefficients  to  account  for  the 
reflections  and  transmissions  undergone  by  the  particiilar  wave  reflection 
In  reaching  (l,j),  is  added  to  the  soam  of  stresses,  S,  already  computed  for 
(l,j),  and  the  computation  proceeds  to  the  next  step;  viz.,  testing  to  see 
whether  the  succeeding  reflection  has  reached  (or  paeBed)(l, J),  If  not^  the 
conqjutation  goes  on  to  the  next  cycle;  it  having  been  found  that  a  given 
reflection  has  not  yet  reached  a  given  point,  there  is  no  reason  to  test 
for  succeeding  waves. 

The  cycle  on  the  left-hand  side  of  Figure  4.7,  and  the  first  three  steps 

on  the  right,  are  common  to  all  four  thickness  ratios  (-^,  the 

computation  goes  throaigh  these  steps  for  all  plates.  In  addition,  there  is  one 

reflection  pecoillar  to  the  plate  (KH=1),  and  two  to  the  plate |at  the 

corresponding  stsige  in  the  lower  plate  (part  )»  there  are  three  for  the 
2 

plate),  so  here  the  program  branches  according  to  the  value  of  KH,  and 
the  indicated  computation  performed. 


Now  the  ret\im  is  made  to  the  main  program,  the  stress  multiplied  by 
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TOl  (the  last  step  in  the  computation^  which^  however#  helng  common  to  © 
and(^^  #  Is  placed  in  the  main  program)  and  printed#  and  the  remainder  of 


the  program  executed. 


4,12.14-  Results  of  Ntnnerlcal  Evfti nation 


The  results  of  this  phase  of  the  program#  as  with  the  single 
plates,  are  presented  in  the  form  of  stress  contours#  figures  4.9  -  4.l4# 
with  the  tensile  areas  shaded  and  the  wave  fronts  indicated  as  in  Figures 
4.4  and  4.5*  The  absence  of  contours  at  the  lower  right  in  several  of  the 
figures  is  due  to  the  uncertainty  resulting  from  the  necessity  of  using  the 
relatively  coarse  mesh  (.lx  .l)  in  order  to  avoid  exceeding  the  capacity  of 
the  con^juter;  however#  the  values  of  the  stresses  in  these  areas  were  small# 
in  general  not  in  excess  of  .05. 

4.13  Qualitative  Discussion  of  Resiilts 

In  order  to  consider  the  relative  merits  of  the  six  laminated  plates 
under  consideration#  they  were  each  divided  into  five  sections  as  shown: 
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The  plates  are  assigned  nvmibers  in  the  same  order  as  that  in  which  they 


appear  In  the  figures: 


1 

.2  Pe, 

.6  Ti 

2 

.4  Pe, 

,k  Ti 

3 

.6  Fe, 

.2  Ti 

k 

.2  Ti, 

.6  Pe 

5 

.U  Tl, 

.4  Fe 

6 

.6  Ti, 

.2  Pe 

Then  for  each  of  the  five  sections,  corresponding  to  the  first  few  wave 
reflections,  the  plates  were  ordered  according  to  the  maximum  tensile  stresses 
fo\xnd  In  the  section.  The  results  are  as  follows: 


Section 

^PP 

I 

1  =  2  =  3>4  =  5  =  6 

(No  tensile  stresses) 

II 

1>2>3>5>6>4 

( None ) 

III 

3>1>2>5>4>6 

( None ) 

IV 

1  >  4  >  5;  None  for  2,3,6 

(None) 

V 

1>2>3>1j.>5>6 

3>2>1>5>I<-  =  6 

In  general,  we  may  make  the  following  ohsei^ations  regarding  the  various 
laminated  plates:  first,  that  those  with  the  Fe  con^jonent  on  top  are  subjected 
to  higher  stresses  j  and  second,  that  those  with  the  thinner  top  plate  are 
also  higher. 

Taking  the  "optimum  plate  to  be  that  with  the  lowest  maximum  tensile 
stress,  then,  we  wo\ild  rate  5^6  first,  then  and  followed  by  ^  and  #3» 
and  finally  which  shows  the  highest  maximum  stresses. 
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4.l4  Implications  on  Ballistic  Resistance  of  Plates 


In  the  previous  paragraphs  of  this  section  the  treuislent  stress  wave 
that  crosses  the  plate  system  has  been  analyzed.  The  two  stress  components, 
radial  (equivalent  to  tangential)  and  vertical  have  been  numerically 
determined. 

In  a  plate  composed  of  a  single  material,  the  two  stress  components  are 
compressive  at  the  wave  front.  The  magnitude  of  these  stresses  Is  quite 
large  and  for  plates  of  thickness  less  than  the  bullet  radius  (h/a<C  l)  the 
vertical  stress  is  of  the  order  of  magnitude  of  the  impact  pressure. 

The  transient  spikes  decay  very  quickly.  The  decay  time  is  of  the  order 


where  a  =  radius  of  bullet 

c  =  longitudinal  sonic  velocity 

The  stress  components  decay  to  a  very  steady  stress  configuration.  The 
vertical  stress  is  compressive  and  the  radial  stress  is  tensile.  These 
stresses  decay  with  depth  much  faster  than  the  spikes. 

Let  us  make  a  quick  estimate  of  the  decay  time  of  a  .1"  steel  plate  for 
h/a  =  .5  (a  ~  0.05  ft).  The  sonic  velocity  for  steel  is  about  17^000  ft/sec. 
Thus 

~  3  X  10'® 

C?) 

This  time  is  very  nearly  the  minimum  to  fracture  observed  by  Zhurkov''  ' . 
Since  this  spike  is  held  for  such  a  short  time,  fracturing  will  not  occur 
unless  the  stress  level  exceeds  the  stress  at  the  wave  front. 
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There  are  three  regions  in  which  the  transient  spike  is  analyzed  for 
fracture : 

1.  As  the  wave  crosses  the  plate,  the  first  time,  the  sheer  stress  or 
the  difference  between  the  compressive  stresses  may  become  large. 

2.  The  reflected  transient  vertical  stress  is  tensile  and  the  resulting 
vertical  stress  may  be  tensile. 

3.  The  reflected  radial  stress  is  tensile  and  the  total  radial  will  be 
tensile . 

In  the  first  wave  the  stresse§*^re  compressive  near  the  wave  front.  The 
fracture  here  will  be  dependent  on  the  resolved  shear  stress.  The  Mohr's  circle 
diagram  of  the  fracture  criterion  is  shown  in  Fig.  ^4-.15.  The  critical  shear 
stress  will  be  taken  as  twice  the  tensile  yield  stress. 

The  shear  stress  in  the  directed  wave  is 


(z/a)  is  the  geometrical  decay  of  the  wave  front  and  is  shown  in  Fig.  4.l6. 

The  maximum  radial  and  vertical  tensile  stresses  in  the  reflected  wave 
are  computed  graphically  by  adding  the  value  of  the  static  stress  at  z/a  to 
the  negative  of  the  stress  at  the  wave  front  at  position  (—  +  z/a.) 

Q, 

The  resulting  maxim\am  tensile  stresses  are  plotted  at  a  function  of  z/a 
parameterized  by  h/a.  The  maximum  radial  stresses  are  shown  in  Fig.  4.17 
and  the  maximum  tensile  vertical  stresses  in  4.l8. 
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T 


Figure  4.15  -  Mohr's  Diagram  Showing  Fractvire  Criterion 

in  Compressive  Wave  Front 
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For  each  plate  thickness  the  vertical  tensile  stress  has  a  maximvim  near 
the  center  of  the  plate.  The  radial  tensile  maximum  is  at  the  surface  near 
the  bullet. 


The  envelope  of  the  maximum  tensile  radial  and  vertical  stresses  are 
plotted  as  a  function  of  the  dimensionless  thickness  h/a  (^ig.  4.19.)  Also 
included  is  twice  the  maximvim  shear  stress  in  the  compressive  wave. 


The  plate  will  break  if  any  one  of  these  stresses  exceeds  the  tensile 
yield  stress.  The  upper  curve  consists  of  the  radial  tensile  curve  for  very 
small  thicknesses  and  the  shear  for  larger  thicknesses.  If  the  plate  material 
is  very  strong  in  compression,  the  upper  curve  for  thicker  plates  will  become 
the  tensile  vertical  stress  envelope. 


Denote  the  upper  bound  by  f(h/a).  The  stress  criteria  for  failure  is 
(now  in  dimensional  form) 


a 

cy 


P  = 


fJ.lC 


V 


S 


a 

^  MiC i-t-MgCg  cy 
50  “  M1C1M2C2  f (h/a) 

The  v^Q  as  a  function  of  plate  thickness  is  shown  in  figure  4.20. 

The  overall  reduction  in  bullet  velocity  may  be  calculated  very  simply 


from  the  force  on  the  bullet . 


where  ra 


V 

F 


mdv 

“cTC 


F 


Mass  of  bullet  =  7I 

L  =  length  of  bullet 


velocity  of  bullet 

force  =  normal  stress  times  area 


TT  as 


TT  a 


Micijugcg 

M1C1+M2C2 


V 


S 
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Figure  4,19  -  Envelope  of  Maximum  Tensile  and  Shear  Stresses 

vs.  Plate  Thickness 
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Plate  Thickness  h/a 


Figure  4,20 


Increase  of  the  v 


50 


with  Increasing  Plate  Thickness 
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The  velocity  is  integrated  from  its  initial  velocity  V  to  its  final  velocity 

s 

Vj^  and  the  force  term  is  integrated  from  zero  time  \intil  the  time  that  the 
plate  breaks  which  is  the  svim  of  the  time  to  fracture  and  the  transient  time 
of  the  wave. 


^trans. 


V. 


Micica 


For  most  cases,  t  is  negligible  compared  with  t 


t  J_  Ail  1 

trans.  =  -r- 


nh 

Cl 


and 


R  _ 


1  - 


M1C2 


/ilC2+M2C2 


h 

L 


where  n 
n 


1  if  failure  by  shear  stress  is  direct  compressive  wave 

2  if  failure  by  reflected  tensile  stresses 


Thus  the  Vj^  versus  curve  is  a  straight  line  which,  if  extended  would 

pass  through  the  origin.  At  the  lower  end  of  this  curve,  however,  the  time 
to  fracture  may  not  be  negligible  and  the  Vp  vs  curve  will  be  slightly 
below  the  predicted  curve. 

The  V  values  and  the  maximum  v„/V  ratios  were  calculated  for  several 
50  K  S 

cases  for  which  experimental  data  is  available.  The  experimental  data  was 

supplied  by  Mr.  A.  L.  Alesi  of  the  Quartermaster  Corps  Research  and  Engineering 

Command,  Natick,  Massachusetts.  The  predicted  v  values  are  given  in  table  4.3* 

50 

The  V-  values  are  evaluated  for  the  compressive  stress-  failure  criteria  in  the 

50 
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direct  wave.  The  values  of  the  yield  stress  used  are  indicative  of  the 
range  of  yield  reported  in  the  literature.  The  plate  descriptions  together 
with  the  experimental  data  is  classified.  The  descriptions  of  the  materials 
and  dimensions  are  given  in  Appendix  B  in  the  same  order  as  Table  4.3.  The 
predicted  values  of  Vgo  are  reasonable  and  reflect  the  \incertalnty  of  the 
measured  values  of  the  yield  strengths. 
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Section  5 


fi 

n 

n 

n 

n 


L. 

L 

U 

II 

1  “ 


MODEL  FOR  SHEAR  FAILURE  OF  THE  PLATE 

The  full  solution  of  the  equations  of  motion  for  impact  involve  the 
cooputation  of  both  elastic  and  inelastic  behavior  in  at  least  two  dimensions 
plus  time.  The  linear  elastic  solution  for  this  problem  is  devilisWy  compli- 
cated.  Only  this  modbl  allows  for  an  analytic  solution;  Since  the  stresses 
are  above  the  static  yield  point  for  most  of  the  deformation  process,  an 
approach  is  desired  that  will  demonstrate  the  effects  of  the  inelastic  behavior. 

To  increase  the  mathematical  description  of  the  physics  involved,  and  yet 
retain  a  fairly  reasonably  sized  problem,  the  geometrical  variations  of  the 
problem  must  be  reduced. 

To  effect  this  geometric  sinplifi cation,  we  note  that  after  the  stress 
wave  has  reflected  at  least  once  the  vertical  variation  in  stresses  is  almost 
zero.  So  the  vertical  variation  of  the  stresses  and  displacements  are 
approximated  by  setting 


Further,  since  most  of  the  motion  of  the  plate  and  projectile  is  vertical, 
we  will  neglect  the  radial  displacement. 

5.1  Historical  Background 

This  type  of  model  has  been  used  by  several  authors  to  analyze  Impact  on 

(i' t) 

thin  plates.  Hopkins  and  Prager  treated  the  deformation  of  the  plate  as 

(9  ) 

pure  bending  and  Craggs  treated  the  problem  as  pure  membrane  stresses 
(  to) 

niyushin'  thought  that  the  deformation  was  pure  shear  and  the  Russian  authors 
Kotchekov^^^^^  Bakhshian^^  )(  3  Sokolov^  used  various  visco-- 

elastic  and  linear  visco- plastic  materials  to  carry  out  explicit  calculations. 
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9*2  MsithwTw^tlceJ.  Model 


Assume  a  cylindrical  biiUet  of  radius  length  L  approaches  a  plate 
of  thickness  h  with  an  incident  velocity  The  densities  of  the  bullet  and 
plate  are  jUi  and  /Za,  respectively.  The  bullet  strikes  the  plate.  Denote  the 
region  of  the  plate  directly  \inder  the  bullet  as  the  plugf  Fig,  (  5.1)  (Assvune 
the  propagation  velocity  for  the  transient  spike  stress  waves  is  infinitely 
fast).  Thus  the  momentum  in  the  bullet  is  transferred  instantaneously  to 
the  bullet  and  plug  system. 

The  mass  of  the  biillet  is 

M  =  TT/ZiR^L  (5.1) 

and  the  mass  of  the  plug  is 

m  5=  TT/iaR^h  (5,2) 

The  moment  Tim  in  the  bullet  before  impact  is 

P  =  MVg  (5.3) 

The  resvilting  velocity  of  the  mass  and  plug  system  is  given  by  the 
conservation  of  linear  momentum 

MVg  =  (M+m)v^  (5.4) 

or 

""o  =  ""s  (5.5) 

Ass-ume  that  the  only  motion  of  the  plate  is  in  the  verticsLl  direction 
and  that  the  only  variation  of  displacement  occurs  in  the  radial  direction 
(i,e.  the  variations  in  stress  across  the  thickness  of  the  plate  are  neglected)^ 
Further,  assume  that  the  membrane  stresses  and  bending,  moments  are  zero.-  ■ 
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This  resultliig  deformation  Is  thus  a  shear  deformation.  A  further  assun^lon 
Is  that  the  plug  Itself  does  not  under  go  ary  f\irther  deformation. 

These  assumptions  have  the  effect  that  only  the  shear  deformation  near 
the  rim  of  the  bullet  will  be  Investigated.  At  the  expense  of  full  geometric 
detail,  the  Influence  of  the  non-elastic  mechanisms  of  deformation  will  be 
brought  to  light.  ! 


Under  these  assumptions,  the  equations  of  motion  In  the  plate  reduce 


to  a  single  equation 

T 

ov  rz  .  rz 

p,2  -  =  -  +  - 

dr  r 


(5.6) 


The  eqviatlon  of  mass  continuity  In  the  plate  is 

dv  _  ^"^rz 
dr  dt 


(5.8) 


One  further  relation  is  needed  to  complete  the  description.  This  equation 
is  the  equation  of  state  if  the  material  is  conservative  or  constitutive 
equation  if  the  material  is  not  conservative.  In  either  case,  this  equation 
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Is  a  relation  between  the  stresses  and  strains 


For  the  plastically  deforming  inedla>  this  constitutive  equation  Is  of 
the  form 


-V  =  -y®  + 

'rz  'rz  'vz 


(5.9) 


where 


0 


•V  =  Total  strain  rate 

•rz 


...e  , 


7 


=  elastic  strain  rate 


rz 

=  plastic  strain  rate 

rz 

For  this  problem,  the  elastic  component  of  strain  rate  will  be  given  by 

the  differentiated  form  of  Hooke's  law. 

T 


= 

'rz 


rz 

G 


(5.10) 


with 


G  =  shear  modvilus 


The  plastic  strain  rate  term  can  be  quite  involved,  depending  on  the 
material  and  boundary  condition  (see  Stecftlon  2  for  details).  For  the  present 
derivatioui  it  is  assumed  that  the  plastic  strain  rate  is  e.  function  of  the 
instantaneous  stress  only.  No  explicit  terms  will  be  inserted  for  this  until 
the  numerical  evaluation. 

Then 


(5.11) 
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Figure  -  Initial  Configuration  of  Plate  and  Bullet 


The  constitutive  equation  is  then 
f 

+  <1>(t  ) 

^rz  or  vz' 


(5.12) 


These  equations  must  he  solved  subject  to  the  following  boundary  and 
initial  conditions. 


Initial  conditions  v  =  0 


V  =  V 


r  >  R 
r  =  R 


t  =  0 


Boundary  conditions  (M+m)  —  =  27TRh  7  r  =  R  t  >  0 


Introduce  the  following  dimensionless  parameters 

'^rz 

§  »  Dimensionless  She9.r  Stress 


7 


7 


rz 


T  =• 


's 

ct 


Dimensionless  Shear  Strain 


Dimensionless  Time 
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Dimensionless  Distance 
Dimensionless  Displacement 


Dimensionless  Velocity 


Sonic  Velocity  In  plate 


In  these  dimensionless  variables  the  equations  reduce  toj 


^  T 
ap  p 


Momentum  conservation 

Continuity 

_  by 
hp  3t 

Constitutive  equation 

^  =  —  +  V'(T) 
hr  St 


with 


»(T)  =  0(-^) 


(5.13) 


(5.li^) 


(5.15) 


‘  (5.16) 


A  word  of  caution  should  he  extended  at  this  point.  The  forin  of  the 
function  f6  which  describes  the  plastic  flow  may  or  may  not  be  amenable  to 
scaling.  The  form  that  shall  be  used  In  later  numerical  calcvilations  is  scalable 
(l.e.  it  can  be  written  in  the  form  j6  =  ^(T)).  However,  many  materials  can  be 
envisoned  where  this  does  not  hold. 

The  boundary  and  initial  conditions  are  ; 

Initial  conditions  v  =  o  P 


: :)’  ■ 
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P  =  1  T  >  0 


Boundary  conditions 


where 


°  St 
m  = 

°  27rR^hG 


Mtm 

55r 


Even  with  simple  forms  of  analytics!,  solutions  are  normal  1  y  Imposslhle, 

so  that  numerical  equations  must  he  used.  However*  a  partial  analytic  solution 
sheds  some  light  on  various  features  of  the  solution. 


The  continuity  eqiiatlon  Is  combined  with  the  constitutive  equation 

~  +  V'(T)  (5,17) 

op  OT 

This  eqviation  with  the  moment\nn  equation  forms  a  system  of  totally  hyperbolic 
Ijartlal  differential  equations.  Thus  this  system  can  be  sln^llfled  by  the 
method  of  Rlemann  characteristics.  The  ulterior  goal  of  this  aralysls  Is  to 
find  a  simplified  solution  for  the  Jump  or  shock  conditions  at  the  wave  front 
that  will  propagate  into  the  medium. 

Add  and  subtract  the  momentum  equation  and  continuity  equations, 

^(T)  (5.18) 

3p  3p  p 

Define  the  characteristics  coordinates  a*  jS 
2a  =  T  +  p 

2^  =  T  -  p 

The  characteristic  jS  =  constant  is  a  line  on  a  sonic  wave  which  propagates 
outward  and  the  line  a  =  constant  Is  a  line  on  a  sonic  wave  which  propagates 
Inward. 
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In  terms  of  these  characteristics  the  equations  above  may  be  written 

dV  dT  T  f  \ 

(T)  along  a=  constant 

dV  dT  T  #  » 

"  p  along  /S=  constant 

The  characteristics  are  the  only  lines  on  which  the  solutions  can  be 
dl scont Inuous . 


Now  consider  the  /3  characteristic  that  passes  through  the  point  />=1,T=0. 
Ahead  of  this  characteristic  the  solution  Is  identically  zero,  ^y  considering 
the  solution  to  any  a-characteristic  that  crosses  this  line  (in  differential 
form) 

dT  -f  dV=  (T-tt(T))  djS  (5.19) 

however  d/3  =  0 


Thus 


-  T‘]  =  -  V^] 


where  the  plus  superscript  refers  to  the  values  on  one  side  of  the  jS-charac- 
teristic  and  the  minus  superscript  refers  to  the  values  on  the  other  side. 
Thus,  at  the  wave  front,  the  jump  in  the  dimensionless  velocity  is  equal  to 
the  negative  of  the  jump  in  the  shear  stress. 


Along  the  /3- character! Stic  the  characteristic  equation  is  siiiq)lified  to 

(5.20) 

or  writing  a  as  a  function  of  p  only 


Thus,  at  the  wave  front  a  simpler  equation  giving  T  as  a  ftmction  of  p  is 
available  together  with  the  relation  T  =  -V 
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5-3  Wumerlcal  Method  of  Solution 


The  eqviatlons  to  be  solved  are 

ay  ^  ,  T 

br  3p  p 

HD 

bp  br 


(5.22) 


(5.23) 


in  the  region 


p  it  1 
T  a  0 


subject  to  the  boundary  and  initial  conditions. 
Initial  conditions  y  =  0  P  ^ 


V  =V, 


Boiindary  conditions 


■&  '• 


p=  1  ,  T  > 


(5.24) 


The  displacement  is  U  and  the  total  strain  and  particle  velocity  are 


defined  by 


V=  M 
br 

7  =  M 


(5.25) 


The  total  plastic  strain  at  any  time  and  position  is  given  by  the  integral 


7„(P,t)  =  /  f{p,T^)  dT» 

^  o 


(5.26) 


In  addition  to  these,  a  simpler  equation  gives  the  stresses  and  particle 
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At 

.-.1 


velocity  at  the  wave  front 

.  l(|+^(T)) 

V  =  -T 


along  p  =  T  +  1 


These  eqiiatlons  were  solved  hy  a  finite  difference  set  of  eq\iatlons. 

The  finite  difference  scheme  is  given  below. 

Form  a  square  mesh  work  in  the  p^r plane.  The  increment  of  mesh  size 
in  T  is  Arand  the  increment  in  p  is  Ap,  as  shown  in  the  figure.  The  index  is 
refers  to  the  increments  in  p  and  the  index  n  refers  to  the  Increments  in  T. 

At  some  time  t“”  the  displacements  and  velocities  at  the  corners  of 
the  mesh  are  known  and  the  stress  and  accumulated  plastic  strain 

at  the  midpoint  between  each  pair  of  mesh  points,  ^  p  i+l/2  known. 

In  addition,  the  wave  front  is  located  at  1  =  n-1  and  the  stress  and  particle 
velocity  at  this  point  are  known.  The  displacement  is  zero  at  the  wave  front. 

First  the  displacements  at  time  are  computed 


for  the  points  1  S  i  S  n-1. 
The  displacement  at  i  =  n  is  zero,  i.e.,  =  0. 


The  total  strain  at  the  midpoint  of  each  cell  is  computed 


7t  n 

^  1-i/s 


^rM-1 


(5.27) 
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The  Increment  In  plaetlc  strain  conqwted  by  a  Rm^e-Kiitta  inte¬ 

gration  assuming  a  constant  tdtal  strain- rate  process  in  each  cell  during 
this  time  step. 

’^p  1+l/Z  “  1+1/2  ^4^  •  (5.28) 


where 


^2  '  "l) 


fi*  =  ^  '■3) 


The  stress  at  each  midpoint  is  fotind  by  the  integrated  Torm.  of  the 
constitutive  equation 


1+1/2  “  ■''t  1+1/2  ■  ■’'p  1+1/2 


(5.29) 


The  accelerations  and  particle  velocities  are  computed  at  each  interior 


point 


^1+1/2 


~  '^1-1/2  1  ^^1+1/2  ^1-1/2^ 

3  'S  p. 


^+1/2  ^  ^-1/2  ^ 


(5.30) 


C5.3I) 


The  stress  euad  velocity  at  the  wave  front  are  computed  by  a  Runge-Kutta 


solution  of 


I;  =  - -5(1  + 


(5.32) 
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The  particle  velocity  at  1  «•  1  Is  solved  by  using  the  boundaxy  condition. 
Interpolation  of  -the  stresses  In  the  first  two  cells  Is  used  to  derl-ve  the 
stress  at  p  s  1. 


IIIq  •^l 


•5  ^2 


tS) 


(5.33) 


All  the  -values  at  have  been  counted.  J^ow-  the  time  Is  Increased  and 
the  cycle  is  repeated. 


The  computation  Is  continued  \intll  one  of  two  things  occurs,  either  the^ 

stress  at  the  surface  goes  below  the  yield  stress,  at  which  time  no  fui^her 

✓ 

/ 

flow  can  occur,  or  the  plastic  strain  in  the  first  cell  exceeds  the  limit, 

at  which  fracture  -will  probably  occur, .  at  which  time  fracture  occurs  and 

the  plug  is  free.  If  the  second  case  occurs  then  the  plate  has  been  perforated. 

The  difference  scheme  was  checked  for  stability  and  convergence  by 

running  trial  calculations  using  different  sizes  of  the  mesh.  The  difference 

k 

scheme  -was  stable,  varying  roiaghly  as  the  order  of  h  . 

5.4  Resvilts 

The  celculations  were  e-valtiated  for  several  cases  of  interest.  The  cal¬ 
culations  are  di-vlded  into  two  categories.  The  first  category  was  the  deter¬ 
mination  of  the  effect  of  the  -various  parameters  on  the  shax>e  and  location 
of  the  Vj^  vs.  Vg  curve.  The  second  category  was  the  e-valtiatlon  of  certain 
specific  plate  and  target  materlsls  for  which  experlmentsl  -values  of  the  Vp-v^ 
were  available.  In  addition  to  these  cases,  one  of  the  cases  -was  e-valuated 
for  the  stresses,  displacements  and  particle  -velocities  e-vezywhere. 
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The  form  of  the  ^aetlc  strain  rate  tenn  In  the  constitutive  equation 
was  taken  as  (see  See.  2) 

6p  -  .  0(t)  (5.34) 

where 

£p^l8  the  maximum  posslhle  plastic  strain  rate 
A  Is  a  constant 

This  form  of  the  plastic  strain  rate  ccmqjonent  Is  easily  rendered  into  the 
dimensionless  fozm  by  appropriate  change  In  A.  Since  the  behavior  of  this 
function  Is  almost  as  sharp  as  a  unit  step  function,  the  scaling  behavior 

■^rz 

is  dependent  on  the  ratio  of  -j— •  For  values  of  this  ntnriber  less  than  0.5, 
the  behavior  Is  almost  entirely  elastic.  For  values  of  this  ratio  In  the 
neighborhood  of  unity  or  above  the  behavior  Is  almost  that  of  a  perfect  fluid. 

t  °  and  A  were  evaluated  for  typical  tltantlum  alloys  by  assianlng  the 
P 

model  derived  In  Section  2, and  uslijg  the  constant  strain  rate  data  for 

(75 

various  tltantlum  alloys'  ' .  ■ 

Table  5*1  gives  the  values  of  A  and  t  The  maximum  critical  strain 

P 

was  ass\jmed  Independent  of  the  rate  and  was  chosen  as  l/2  of  the  critical 

(35) 

elongation  in  tension.  The  most  probable  value  In  the  literature  Is  also 
shown  In  the  table. 


The  residual  velocity  v„  versus  the  striking  velocity  v„  was  evaluated  for 

X\  8 

two  typical,  tltantlum  alloys  which  show  the  Influence  of  one  of  the  parameters. 
The  two  alloys  were  Ti  2.5  Sn  5A1  and  Ti  4  Mn  4a1  which  very  nearly  have  the  same 


static  conpresBlve  yield  strength,  the  same  value  of  A,  and  the  same  value  of  e 
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Table  5-2 

MATERIAL  PROPERTIES  OF  PLATE  ARP  PROJECTILE  FOR  DETAILED  CALCPLATIOWS 


Property 

Pirojectile 

Plate 

Material 

43160  Steel 

T1  2,5  Sn  5  A1 

Pertinent 

Dimensions 

R  =  0.22" 

L  =  0.25" 

h  =  0.100" 

Density 

0,283  Ibs/in 

0.163  Ibs/in^ 

istatic 

compressive  yield  strength 

150','000  psi 

120,00  psi 

A  coefficient  in  exponent 
of  strain  rate  term 

8,0  x  10^  psi 

0 

e  max  plastic  strain  rate 

r  ■ 

5  X  10^/sec 

max  critical  elongation 

0.15 
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The  major  difference  Is  that  the  critical  i>emanent  strain  at  fracture  In 

one,Tl  2.5  Sn  5Alyls  nearly  double  that  of  the  other  alloy.  The  Vp  vs.  v^  curves 

that  were  calculated  for  these  two  cQ-loys  are  shown  In  figures' 2  -  5*9« 

The  v^Q  point  or  the  velocity  of  incident  perforation  is  taken  as  the  Intercept 

of  the  v_-v  curve  with  the  v„  axis.  This  v^^  is  plotted  versus  m  for  the 
R  s  s'  50 

two  alloys.  (Fig.  5-10)» 

A  typical  steel  was  the  other  case  evaluated.  The  same  yield  strength 
as  the  tltantl'um  alloys  was  used  but  the  maximum  strain  rate  and  maximum 
critical  elongation  were  taken  as  twice  that  of  the  Ti  2.5  Sn  5A1.  No  good 
strain  rate  data  Is  available  on  steels.  This  data  was  chosen  since  It 
represents  a  much  more  ductile  material,  The  v„  vs.  y  curves  are  shown  In 
figures  5*11“5.1^.‘  ‘  '  ■  . 

One  case  was  evaluated  for  all  stresses  and  particle  velocities  as  a 
f\inction  of  position  and  time  for  several  values  of  the  striking  velocity. 

The  case  chosen  for  this  evaluation  was  a  Ti  2.5  Sn  5A1. plate  thickness  0.10  inches 
being  struck  by  a  43^0  steel  cylindrical  projectile  0.22"  diameter,  O.25  " 
length.  Material  properties  of  the  plate  are  given  in  Table  5j2.  The 
vertical  particle  velocity  and  the  shear  stress  are  shown  as  a  function  of 
position  and  time  for  five  impacting  velocities,  ranging  from  below  the  elastic 
limit  to  25  times  higher,  '  (Fig* 9 *15-24).  Where  the  diagrEims  are  imcomplete,this 
means  the  plate  has  been  perforated  and  the  stresses  have  gone  to  zero. 

The  permanent  deformation  of  the  plate  as  shown  for  these  five  cases 
in  Figure  5.25. 

Four  specific  cases  of  the  v„  vs.  v  cxirves  were  evalviated  to  correspond 

K  8 
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Figure  5^6 

Vp  vs  Vg  For  a  T1  Mn  4  A1  Target 

with  m  =  1 
o 
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Figure  5.7 

V„  vs  V  For  a  T1  4  Ml  4  A1  Target 

A  8 

with  m  =  2 
o 


MRD  DIVISION 

OENERAI.  AMERICAN  TRANSPORTATION  CORPORATION 


CQ 


CO 


MRD  DIVISION 

OENERAI.  AMERICAN  TRANSPORTATION  CORPORATION 


109 


0-j - 1 - 1 - 1 - 1 - 1 - 1 - 

0  123  ^5  6  **0 

Figure  5»10  -  vs.  m  for  Two 

50  o 

Titanium  Alloy  Targets 
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2.3  2.5 


Figure  5 *25  -  Permanent  Vertical  Displacement  as  a  Function  of  Radius  and  Input  Velocity 

MRD  DIVISION 

OBNERAI.  AMERICAN  TRANSPORTATION  CORPORATION 


125 


with  experlmentEil  measiirementB  of  vs.  v^.  This  experimental  data  was 
supplied  "by  Mr,  A.L.  Alesl  of  the  Quartermaster  Corps  Research  and  Engineering 
Command.  The  description  of  the  plate  and  projectile  materleO-s  In  conjunction 
with  the  experiments  Is  classified.  These  descriptions  of  the  plate  and 
projectile  materials  are  given  In  a  classified  Appendix  B. 

The  v„  vs,  V  curves  together  with  the  experimentally  determined  points 

X\  8 

are  shown  In  Flg.5‘-26-28.  The  agreement  "between  the  theoretically  predicted 
results  and  the  experimental  resxxlts  Is  considered  good. 

5.  J  Summary 

fhls  section  has  investigated  a  simple  geometrical  model  of  penetration 
using  a  physically  nonlinear  model  of  nonelastic  behavior, 

The  stress-position-time  pattern  is  altered  if  there  is  plastic  flow. 

There  are  three  basic  parameters  in  the  nonelastic  model  that  Influence 
the  Vp  versus  v^  curves,—  he  exponent  coefficient  A (which  is  directly  related 
to  the  static  yield  strength),  the  maximum  flow  rate  and  the  critical  strain 
for  fracture,  In  addition,  there  is  the  geometrical  factor  of  thickness  or 
total  mass  of  the  plate  under  the  bullet. 


For  all  the  material  properties  chosen,  if  these  are  held  constant  and 


only  the  thickness  of  the  plate  is  increased  the  v^^  decreases  almost  linearly 

,  —  Jit-m 
“o  = 


If  the  flow  rate  of  the  id.ate  is  held  constant, 
increasing  critical  elongations. 


the  V 


50 


Increases  with 
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2000  3000  hOOO 


Since  in  most  materials  the  floir  rate  emd  the  critical  elongation 

Increase  together  (this  coniblnatlon  Is  most  likely )>  one  set  of  v„  versus 

V  curves  was  evaliiated  Increasing  both  by  the  same  percentage.  In  this 
s 

case  the  v^q  or  resistance  to  penetration  was  significantly  reduced.  The 
stresses  In  the  region  of  a  large  amount  of  plastic  flow  are  reduced  slower 
than  if  there  is  little  flow.  Thus  the  plastic  strain  at  that  point  will 
increase  faster,  thus  tending  to  break  the  specimens. 

The  Vp  vs  Vg  curves  for  high  impact  velocities  are  straight  which,  if 
extended,  would  pass  throiigh  the  origin.  This  reflects  the  momentTom  transfer 
to  the  plate  material  proportional  to  the  Impact  velocity.  This  momentum 
transfer  occurs  over  a  short  duration  of  time  which  is  the  time  to  fracture 
at  very  high  stress.  At  lower  impact  velocities  the  stress  in  the  plate  is 
lower,  the  time  to  fracture  becomes  loiiger,  more  momentum  is  transferred  in 
the  plate  dtiring  this  longer  interval  and  the  vs  Vg  curve  has  a  greater 
shape  and,  if  extended,  will  intersect  the  Vg  axis  to  the  right  of  zero.  The 
ciorvatvire  in  this  bend  is  not  great,  but  the  intercept  is  shifted  by  anywhere 
for  50  fps  to  500  fps. 
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Section  6 

COMBIHATION  OF  RESULTS 

In  the  previous  two  sections^  two  models  describing  two  critical  phenomena 
In  the  perforation  of  metallic  plates  and  metsilllc  laminates  have  been  discussed. 
The  two  critical  areas  are: 

1.  Transient  wave  propagation 

2»  Long  time  (con^jared  to  the  transient)  plastic  flow 

The  Vp  vs  Vg  ctirves  and  the  v^q  I)olntB  were  predicted  In  each  case.  In 
the  actual  problem  both  problems  must  be  discussed  simultaneously,  both  mechanisms 
of  perforation  of  the  plate  must  be  considered  In  determining  the  resistance 
of  the  plate  or  laminate  to  ballistic  Impact, 

The  sln^jlest  way  of  showing  the  combined  behavior  Is  to  assiune  that  each 
mechanism  Is  entirely  Independent  of  the  other  and  then  to  superpose,  the  v^ 
vs  Vg  curves.  The  lower  v^^  and  the  higher  v^  over  Vg  ratios  must  be  used. 

This  combination  of  curves  Is  shown  schematically  In  fig.  6.1  for  a 
single  plate.  From  the  res\ilts  of  sections  k  and  5  the  transient  wave  Is 
normally  the  dominant  mode  of  failure. 

In  the  laminated  plates,  the  transient  wave  may  be  reduced  In  the  rear 
plate  so  that  the  critical  model  for  failure  will  be  the  shear  model.  This 
Is  shown  schematically  In  fig.  6.2, 
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Figure  6,2 


-  Vp  vs  V„  for  Laminated  Plate  (Schematic) 
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Section  7 


COHCLUSIOHS  AKD  RECX)MME!WDATIONS 
7.1  Survey  of  Results 

The  following  sequence  of  events  takes  place  when  a  biillet  or  projectile 
strikes  a  plate  or  combination  of  plates.  When  the  bullet  strikes  the  plate 
a  stress  is  immediately  established  at  the  interface  between  the  two.  A 
stress  wave  is  propagated  into  the  plate  and  a  stress  wave  is  propagated 
back  into  the  biollet. 

The  stress  wave  that  propagates  into  the  plate  has  very  high  compressive 
stresses  at  the  wave  front.  Both  the  radial  and  vertical  stress  are  compressive. 
In  a  very  short  time  after  the  wave  front  has  passed,  these  stresses  are 
altered.  The  vertical  stress  is  still  compressive  but  reduced  in  magnitude. 

The  radial  stress,  however,  changes  size  and  becomes  tensile. 

The  stress  wave  in  the  plate  traverses  the  plate  and  is  reflected  from 
the  back  surface.  Both  of  the  reflected  stresses  at  the  wave  front  are  tensile. 
These  stresses  add  to  the  existing  stress  from  the  tail  of  the  downward  wave. 

In  most  cases  the  vertical  stress  is  very  nearly  zero  since  a  compressive 
component  and  a  tensile  component  are  being  added.  However,  the  end  result 
is  a  tensile  vertical  stress  which  results  in  spalling. 

On  the  other  hand  the  radial  stress  in  the  tail  of  the  downward  wave  is 
tensile  and  the  reflected  radial  stress  at  the  wave  front  is  also  tensile. 

When  these  two  add,  the  resulting  radial  stress  is  large  tensile. 

Due  to  this  transient  behavior  of  the  propeigatlng  stresses  the  plate 
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cam  be  broken  in  three  vays: 

1.  By  speJ-llng^  due  to  reflected  tensile  vertical  stress. 

2.  By  radial  and  tangential  fractvjre  due  to  the  large  tensile  radial 
stress  of  the  Incident  and  reflected  radial  stress. 

3.  By  compressive  shear  fracture  as  the  first  wave  crosses  the  plate. 

If  the  plate  survives  the  first  wave,  the  problem  of  perforation  is 
not  over.  The  bullet  has  been  reduced  somewhat  In  velocity  but  it  is  still 
moving.  Now  this  process  of  fracture  is  akin  to  extrusion  or  punching^.' 

The  dominant  mechanism  is  now  shear  failure  at  the  rim  of  the  projectile. 

The  laminated  plates  after  the  transient  stress  wave  characteristics 
by  internal  transmission  and  reflection  of  stresses.  The  proper  choice  of 
lamlnatloii  materials  will  reduce  the  critical  stresses  in  the  back  plate. 

The  effects  of  lamination  may  be  summarized  by  considering  the  transmission 
and  reflection  properties  of  the  lameller.  There  are  two  distinct  cases, 
based  on  the  order  of  the  plates.  The  first  case  is  where  the  upper  plate 
has  a  high  mechanical  impedance  and  the  lower  plate  has  a  low  mechanical 
impedance.  The  second  case  is  where  the  plate  order  is  reversed. 

In  either  case,  the  stress  wave  transmitted  to  the  back  plate  is  lower. 

The  reflected  wave  is  tensile  in  the  first  case  and  compressive  in  the  s^cjQj^d'. 
The  second  reflection  of  this  wave  Is  off  the  original  surface  and  is 
con5>ressive  and  tensile,  respectively.  In  either  case  large  tensile  stresses 
develop  in  the  upper  plate.  The  tensile  fracture  is  delayed  a  little  in  the 
second  case. 

The  effect  of  lamination  may  be  summed  up  to  reduce  the  stresses  in  the 
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'back  plate,  the  stresses  In  the  upper  plate  will  be  Increased.  The  upper 
plate  will  fracture  then  at  a  lower  Impact  velocity.  The  lower  plate's 
resistance  to  fracture  by  the  stress  spike  will  be  Increased. 

However,  the  back  plate  must  now  resist  fracture  by  shear  deformation. 

The  projectile  is  moving  at  a  lower  velocity.  The  analysis  of  Section  5  now 
applies  to  the  rear  plate  with  a  lower  incident  velocity. 

7.2  Design  Implications 

The  objective  of  this  study  progrsun  is  to  make  a  theoretical  determination 
of  the  phenomena  that  occur  when  a  projectile  strikes  a  piece  of  armor. 

Through  a  knowledge  of  such  phenomena,  the  factors  which  deter  the  penetration 
of  the  armor  can  be  optimized  and  more  effective  armor  designs  can  be  realized. 
Although  this  study  does  not  complete  the  total  theory  required  to  optimize 
armor  design,  we  can  already  draw  certain  conclusions  which  will  help  the 
armor  designer. 

Results  of  the  study  up  to  this  point  indicate  that  there  are  two 
effects  which  must  be  eliminated  or  significantly  reduced  in  order  to  enhance 
the  ability  of  personnel  armor  to  do  its  Job.  These  effects  are  due  to 
transient  stress  spikes  and  long-term  shear  stress. 

The  effect  of  a  stress- spike  perpendicular  to  the  plane  of  the  plate  is 
spallation.  The  effect  of  a  stress- spike  in  the  plane  of  the  plate  (radial 
stress)  is  tearing  of  the  plate.  These  effects  may  be  alleviated  by  using 
laminated  plates.  The  degree  of  alleviation  will  depend  on  the  materials 
Euad  geometries  that  are  feasible  to  use  within  the  limitations  imposed 
(economics,  weight,  availability,  bulk,  etc,). 
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£^ren  If  the  effects  of  transient  stress  spikes  are  significantly  reduced# 
however#  the  long-time  shear  stress  Is  a  factor  still  to  he  considered.  It 
Is  well  known  that  the  effects  of  long-time  shear  stress  can  he  overcome  hy 
Increasing  the  thickness  of  the  armor.  Since  Increasing  the  thickness  of 
present  armor  designs  would  also  Increase  their  weight#  this  approach  Is  not 
desirable. 

Thus#  the  In^rovement  of  armor  Is  somewhat  dependent  on  a  more  adequate 
use  of  materials  research.  With  the  results  of  the  analytic  program  to  date 
and  with  anticipated  results  of  future  analytic  work#  a  straight-forward 
and  direct  approach  should  he  possible, 

7.2.1  Materials 

Conclusions  reached  from  the  analytic  studies  to  date  should 
he  helpful  in  providing  guidance  in  the  search  for  new  materials  to  use  for 
protective  armor.  This  guidance  may  he  used  in  two  ways. 

The  capabilities  of  new  materials  as  protective  armor  can  he  assessed 
hy  making  a  theoretical  calc\ilation  of  This  would  eliminate  the 

need  for  experimenting  with  materials  which  offer  no  hope  for  improvement. 

In  the  development  of  new  material  concepts#  the  theoretical  conclusions 
may  he  used  to  set  the  direction  of  research.  The  qualities  which  enhance 
a  material's  capability  to  resist  penetration  become  the  objects  of  iii5)rovement . 
Conq)Osite  materials  may  he  developed  which  combine  the  best  qualities  of  two 
or  more  materials  to  produce  a  generally  better  armor  material.  For  exaaj^le# 
it  Is  known  that  many  ceramic  materials  possess  unusually  high  compressive 
strength;  however#  they  are  weak  in  tension.  It  woiold  seem  that  a  fruitful 
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avenue  of  investigation  would  lie  in  developing  a  con^slte  material  which 
would  make  use  of  the  desirable  properties  of  the  ceramic  by  combining  it 
with  such  a  material  as  woiild  overcome  its  weaknesses.  Since  several  material 
properties  are  involved  in  resistance  to  penetration,  the  problem  will  be  one 
of  optimization.  Such  a  composite  might  take  the  form  of  a  ceramic  "slab" 
reinforced  by  high  tensile  strength  wire  in  a  suitable  configuration. 

7.2.2  Geometry 

In  addition  to  the  Improvements  that  may  be  obtained  solely 
through  materials  research,  there  are  indications  that  improvements  are  possible 
through  geometry  as  well.  Atkins^ has  reported  an  experiment  wherein 
three  separated  thin  plates  were  more  effective  against  penetration  by  a 
hypervelocity  fragment  impact  than  a  single  plate  of  more  than  5  times  the 
combined  thickness  of  the  spaced  plates.  A  photograph  of  the  results  of  his 
experiment  shows  that  the  projectile  neatly  punched  through  the  first  plate 
encountered,  ripped  a  large  jagged  hole  and  severly  bent  the  second  plate, 
and  bent,  but  did  not  penetrate,  the  third  plate.  It  would  seem  from  these 
results  that  the  first  plate  deforms  the  projectile  so  that  it  presents  a 
rather  large  blunt  surface  to  the  second  plate.  The  second  plate  absorbs 
much  of  the  energy  from  the  projectile,  thus  greatly  reducing  its  velocity. 

The  third  plate  was  then  able  to  resist  penetration. 

Atkins'  experiment  used  a  separation  of  4  inches  between  plates.  It 
is  felt  that  this  much  separation  is  unnecessary.  The  only  spacing  that  is 
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necesseoy  'between  the  first  and  second  x)lates  Is  that  required  to  allow  the 
projectile  to  deform  before  reaching  the  second  plate.  The  spacing  between 
the  second  and  third  plates  shoiild  be  sufficient  to  allow  them  to  act  on  the 
plate  independently. 

Each  of  the  three  plates  should  be  designed  of  such  material  as  to  maximize 
the  effect  it  is  to  produce.  The  first  plate  should  be  of  a  material  that 
will  cause  greatest  deformation  of  the  projectile.  The  second  plate  shotild 
be  designed  for  maximum  energy-absorbing  damage  to  itself  in  order  to  reduce 
the  velocity  of  the  deformed  projectile  as  much  as  possible.  The  third  plate 
should  be  designed  to  withstand  the  penetration  of  the  deformed,  reduced- 
velocity  projectile.  In  application,  the  spacing  between  the  plates  could 
be  maintained  by  a  suitable  light-weight  material  (e.g.  cotton  batting). 

7.2,3  Local  Effects 

It  should  be  remembered,  in  any  personnel  armor  design,  that  the 
problem  of  penetration  is  local.  The  static  configuration  considered  as  a 
structure  is  not  of  significance.  Only  those  effects  which  apply  to  the 
immediate  area  of  Impact  shoiild  be  considered. 
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APPEHUDIX  A 


FORMAL  INTEGRAL  SOLUTION  FOR  STRESSES 
IN  A  LAMINATED  RLATE 


■by 

J.  D,  Stein  and  P.  G.  Boekhoff 
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A.  I  iDtroductlon 


The  problem  of  finding  formal  expressions  for  the  stresses  In  a  laminated 
plate  will  be  investigated.  It  Is  assumed  herein  that  both  the  material 
properties  of  the  plate  and  the  otherwise  arbitrary  load  on  the  plate  are 
auclally  symmetric.  The  plate  will  be  loaded  dynamically  with  an  arbitrary 
pressure  p(r,t),  and  a  method  for  finding  the  stresses  in  the  plate  in  terms 
of  the  loading  function  and  the  plate  characteristics  will  be  developed. 


I 


First,  the  potential  fiinctions  0  and  6  are  defined  and  the  partial- 
differential  potential  equations  are  presented.  To  these  equations  and  the 
appropriate  boundary  conditions,  Hankel-Laplace  transforms  are  applied.  A 
set  of  equations  in  the  transform  space  is  developed  which,  when  solved  by 
determinants,  yields  expressions  for  the  transformed  potential  functions 
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and '7.  Inversion  foxmilae  are  given  by  which  0  and  0  can  be  found,  and 
the  stresses  finally  appear  as  explicit  functions  of  the  potentials  0  and  6  . 


A. II  Eqviatlons  of  Motion 

Under  the  assumption  of  axial  symmetry,  the  equations  of  motion  can  be 
written  as 

br  a  - 

p — E _ 5:  +  _i2  +  _j: — S 

Sr  Sz  r 


St 

+ 

=  . 

Sr 

Sz 
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St 

rz 

ba 
+  _ z 

Sr 

Sz 

and  u 

are  the 

radial 

and  a^,  Oqj  a^,andT^^  are  the  radial,  tangential,  axial,  and  shear  stresses 
respectively,  and  p  Is  the  plate  density.  The  stresses  are  then  given  by 


r 

Sz 
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Su 

Su 

Su 

t 
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T.  =  p  ( - -  +  - -) 

br  ^ 


where  X  and  p  are  the  Lame  numbers  of  the  plate  material. 
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A .  HI  Boundary  Conditions 


At  the  surfaces  of  the  laminated  plate,  the  boundary  conditions  may  be 
slii5>ly  formulated  and  are  not  dependent  on  the  characteristics  of  the  plate; 
at  the  Interface,  however,  some  assumption  must  be  made  regarding  the  charac¬ 
teristics  of  the  bond.  At  one  extreme,  the  bond  may  be  considered  Infinitely 


strong,  in  which 

u  (1) 

case  we  have 

=  u  (2) 

z  = 

h 

r 

r 

1 

u  (1) 

II 

z  = 

h 

z 

z 

1 

du 

r 

II 

ro 

z  = 

h 

dr 

dr 

1 

au  ^2) 

_  _ z 

z  = 

h 

dr 

dr 

1 

where  the  superscripts  and  denote  the  top  and  bottom  plates  respectively. 
At  the  other  extreme,  taking  the  bond  strength  to  be  zero  gives 


<1)  =  T  (2)  =  0 
rz 


rz 

a  g  0 


z  =  h 


z  =  h 


(J  (2)  <  Q 
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(2)  =  0 


u  (1)  =  u  (2) 
z  z 


z  =  h  ,  C  'Cr 

z  =  h  j  0 

X  ^ 

z  =  h  ,  ^  0 

1  t 
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z  z 


z  =  h  ,  <  0 

1  z 


dr  ® 
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We  shall  assvune  here  a  bond  such  that  the  vertical  stresses  and  displace¬ 
ments  are  eqtial  across  the  Interface  and  the  rz- shear  stress  Is  zero  at  each 
side  of  the  Interface.  Then  If  p(r,t)  Is  the  pressure  on  the  top  surface^ 
the  boundary  conditions  are: 


=  -p(r,t) 

a  (1)  =  p  (2) 
z  z 

a  =  0 
z 


z  =  0 

z  =  h 
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z  =  h 
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rz 


=  0 


z  = 
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rz 


(2) 
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z  = 


h 
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(1) 


z  = 


h 
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A.  IV  The  Potential  Equations 


The  partial- differential  equations  which  must  be  solved  subject  to  the 
above  boundary  conditions  are 


1  =  1,2 


(A.l) 
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7^  e,  -  li  =  JL- 
^  at® 


\jiiere 


7-  + 

dr®  r  ar  az® 
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2  ^1  •*•  ^^3 


1  =  1,2 


(A.2) 
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0^  and  0^  are  defined  by 


u. 
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(1)  ,  .  l!i 

ar  az 

(1)  =  +  ^  +  ^ 

’  az  r  ar 


(A. 3) 
(A.4) 


and  the  stresses  are  given  in  terms  of  0^  and  0^  by 
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(A.8) 
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If  ve  make  the  substitutions 


T  =  a  t 

1 

\  =  1 

a 

-y  =  — 

a 


a 

a 


V  = 
2 


=  f: 


equations  (A.l)  and  (A. 2)  become 

2  2 
'^1  ^^2 


e. 

v^e.  — -  =  V 


2  _ i 

br^ 


(A. 9) 

(A. 10) 


A.V  Transformation  of  the  Differential  Equations 

We  shall  denote  by  0.  the  Hankel  (order  zero)  -Laplace  transform  of  0^  f 

X  ^  O  ^ 

and  by  ^  the  Hankel  (order  one)  -Laplace  transform  of  .  Multiplying 
equation  (A. 9)  by  r  J^(4r)e”^”^  and  integrating  over  t  and  r  from  0  to»»  we 


obtain : 


-4^ 


dz^ 


‘^1,0 


(A. 11) 
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f! 


Multiplying  (a. 10)  by  r  J  (4r)e"  and  integrating: 

1 


Ij  1  2  aZ  ~S 

r^=^i  ®  ^1,1 


(A. 12) 


Now  we  let 


mi^  =  e  +  7i^ 
so  that  (A.II)  sund  (A. 12)  become 


-  "‘I'^l.o  = 


(A. 13) 


^  =  ° 


(a.i4) 


and  we  obtain  the  following  equations  in  the  transform  space: 

_  m  z  -m  z 

0  =Ae^+Ae^ 

1  j  O  1  2 


(A. 15) 


__  m  z  -m  z 

0  -Ae^-t-Ae^ 

2  >  O  3  4 


(a.i6) 


n  z  -n  z 

U'  =Ae^+Ae^ 

a»  1  5  6 


(A. 17) 


n  z  -n  z 

S’  =Ae^+Ae^ 

2^1  7  8 


(A.18) 


A. VI  Boundary  Conditions  in  the  Transform  Space 

To  evaliiate  the  constants  in  these  last  four  equations,  the  boundary 

conditions  must  be  transformed.  Since  these  involve  only  and 

z  rz 
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we  can  confine  ourselves  to  equations  (A.U),  (A. 7)*  and  (A. 8). 
z 

Multiplying  (a. 4)  by  r  J^(4r)e  and  Integrating  over  r  and  T  from  0  to 
we  have: 


u(i)  =  ^ 

z,o  oz  ^  1,1 


(A.19) 


Multiplying  (A.?)  hy  r  J^(^r)e~  and  integrating 

+  2m^)  7i^  +  2/i^e 


(A. 20) 


Multiplying  (A. 8)  by  r  J  (4r)e”  and  integrating 


1  =  (2^ 
rz,l 


(A. 21) 


Now  let 

Pi  =  (X^  +  2Mj_)  7^^ 

so  that  equations  (A. 19)  -  (A.2l)  become 
^z,o  "  dz  ^  ^1,1 


(A. 22) 


+  q, 

z,o  ^1^1,0  ^1  dz 

/ .  \  d!^ .  -  d^j 

=  -q,  -  fi,  - ^ 


(A.23) 


(A, 24) 


Denoting  by  p^  s)  the  Hankel  (order  zero)  -Laplace  transform  of 
Po(r,T),  we  have  for  the  transformed  bo\indary  conditions : 
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a(l)  =  a(2) 
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Substitution  of  equations  (A. 13)  -  (A.l8)  and  (A. 22)  -  (A. 24)  into  the 


above  conditions  gives  us  the  following  system: 
pA  +pA  +qnA  -qnA=-p 

1  1  1  2  115  1  1  6  0^0 

m  h  -m  h  n  h  -n  h 

pe^^A  +  pe  ^^'A  +  qne^^A  -  qne  ^*-A  = 


(A.25) 


1  1 


2  11 


5  11 
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pe^^A  +  pe  ^^A  +  qne^^A  -  qne  ^^A 
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(A. 26) 


m  h  -m  h  n  h  -n  h 
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(A. 27) 


qmA  -  qmA  +  fi(n^  +  ^^)  A  +  /j.  (n  ^  A  =  0  (A. 28) 

111  112  l'  1  '  5  l'  1  6 


m  h  -m  h  n  h  .  „ 

q  m  e  ^  ^A  -  q  m  e  *•  *-A  +  p  (n  ^  +  4^)e  ^  *-A  +  /i  (n  ^  +  4^) 


2  11 


11  1  11  211  511 

-n  h 

e  ^  *-A^  =  0  (A.29) 
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(A.30) 
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A. VII  Solution  of  the  Equations  in  the  Transform  Space 


(A. 32) 


Let  us  introduce  the  following  notation; 


a  =  p 

11  1 


in  h 

^  2  2 


-m  h 

.  2  2 


m  h 

^  1  1 


q  n  e 
2  2 


n  h 
2  2 


a  = 

22 


-m  h 

.  1  1 


q  n  e 
2  2 


-n  h 

2  2 


a  = 

23 


m  h 

2  1 


a  =  q  m 

41  11 


a  = 

24 


-m  h 

.  2  1 


q  n  e 

1  1 


n  h 

1  1 


q  in  e 

1  1 


-n  h 


q  n  e 

1  1 


-m  h 


q  m  e 

1  1 


q  n  e 
2  2 


n  h 

2  1 


n  h 

|Li  (n  *■  *■ 


1  1 


q  n  e 
2  2 


-n  h 

2  1 


-n  h 

/i  (n  ^+e^)e  1  1 
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a  s 

63 


q  ra  e 
2  2 


m  h 
2  1 


L  =  m  e 

sx  1 


m  h 

1  X 


a  = 

64 


q  m  e 
2  2 


-m  h 

2  X 


a  =  m  e 

82  X 


-m  h 

X  X 


67  2  2 


n  h 

/i  (n  2  ^ 


a  =  m  e 

83  2 


m  h 

2  X 


68  2  2 


-n  h 

ju  (n  ®  ^ 


a  =  m  e 

84  2 


-m  h 

2  X 


a  = 

73 


q  m  e 
2  2 


m  h 
2  2 


85 


n  h. 


L  =  q  m  e 

74  2  2 


-m  h 
2  2 


77  2  2 


n  h 

fi  (n  2  2 


a  = 
86 


a  = 

87 


-n  h 

1  *■ 

n  h 

4e  2  1 


a  = 


78  .22 


-n  h 

fx  (n  ®  ^ 


a  =4® 


-n  h 

2  X 


88 


Now  let  r  denote  the  8x8  matrix. 


a 

a 

0 

0 

a 

-a 

0 

0 

11 

11 

15 

15 

a 

a 

-a 

-a 

a 

-a 

-a 

a 

21 

22 

23 

24 

25 

26 

27 

28 

0 

0 

a 

a 

0 

0 

a 

-a 

33 

34 

37 

38 

a 

-a 

0 

0 

a 

a 

0 

0 

41 

41 

45 

45 

a 

-a 

0 

0 

a 

a 

0 

0 

51 

52 

55 

56 

0 

0 

a 

-a 

0 

0 

a 

a 

63 

64 

67 

68 

0 

0 

a 

-a 

0 

0 

a 

a 

73 

74 

77 

78 

a 

-a 

-a 

a 

a 

a 

-a 

-a 

81 

82 

83 

84 

85 

86 

87 

88 
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and  A  the  coltunn  matrix 


0 

0 

0 

0 

0 

0 

0 

and  r .  the  8x8  matrix  formed  by  substituting  A  for  the  Jth  coliimn  of 

U 

Then  If  1  1  FI  [and  [  [Tj!  |  represent  the  determinants  of  F  andFj  respectively 
we  have 

Aj  =  llr^  ll/llrll  J  = 

The  determinants  are  evaluated  by  standard  matrix  manipulations;  the 
following  sample  evaluation  will  give  an  indication  of  the  method  and  of  the 

form  and  the  degree  of  complexity  of  the  coefficients  A. : 

J 


a 

11 

a 

11 

a 

15 

-a 

15 

0 

0 

0 

0 

a 

41 

-a 

41 

a 

45 

a 

45 

0 

0 

0 

0 

a 

51 

-a 

52 

a 

55 

a 

56 

0 

0 

0 

0 

a 

21 

a 

22 

a 

29 

-a 

26 

-a 

23 

-a 

24 

-a 

27 

a 

28 

a 

81 

-a 

82 

a 

85 

a 

86 

-a 

83 

a 

84 

-a 

87 

-a 

88 

0 

0 

0 

0 

a 

33 

a 

34 

a 

37 

-a 

38 

0 

0 

0 

0 

a 

T3 

-a 

74 

a 

77 

a 

78 

0 

0 

0 

0 

a 

63 

-a 

64 

a 

67 

a 

68 
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A-12 


! 

I 

I 

n 

« * 

fi 

[. 

! 

1 


I 

0 


W 


^11  ° 


b  -2a 


K>  -a  -a 

■>  ai 


K)  a;  - 

aa 


K)  -a  -a 

sa  ai 


k)  0 


a  [a  1 

ai  [  aa  a^^  „J 

“] 


0 

2a 


0 

0 


0 

0 


0 

0 


0 

0 


a  +a  ,  0 

as  b4 


..tBIg  I  a  -a 
^11 


p  -  I  a  +a 
^43^  ail  88 


0  0 

0  0 


r  -1 

a; 

» •  “1 

['^**"^68  »“ 

-a 

L  “*  ®48  “J 

“J 

'a  ■ 

L  ®*  ^68  **J 

'  a  ■  1 

r  a  1 

r  *7’*  1 

r  a  -fy^a 

a  - 

a  -K-feTe  1 

L  ••^68»*J 

L  **  *68  •*] 

L  *68 

r  A  1 

r  1 

I  a 

-a 

L  *68 

L  *68  ■^*J 

.a  ' 

0 

0 

0 

0 

0 

0 

0 

0 

-2a 

41 

0 

0 

0 

0 

0 

0 

0 

0 

Id 

88 

^6 

0 

0 

0 

0 

0 

0 

tj 

88 

^6 

b 

as 

b 

84 

0 

0 

0 

0 

b 

SB 

^6 

b 

88 

b 

84 

0 

0 

0 

0 

0 

0 

b 

88 

b 

84 

0 

0 

0 

0 

0 

0 

0 

0 

b 

T7 

0 

0 

0 

0 

0 

0 

0 

0 

*68 
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A-iS 


where 


a  ra  -  (a  /a  )a  i 

b  ,,  a  -  a  +  \-^ - (a  -  a  ) 

25  25  a^^  21  “41  J 


b  =  a  -  a  +  ^  (a  —  a  ) 

26  25  26  a  '  22  21 

4  X 


a  ra  -  (a  /a  )a  T  a 

^  ^  K':  -  <y<>Zr^^  * 

a  ra  -  (a  /a  )a  i 

a  28  [a  (a  /a  }a  J  27 

68  L  77  67  68  78-1 

a  ra  -  (a  /a  )a  1 

63  „  l  73  63'.  68'  78|/« 

^  ~  ^38  a^ — -  Ca  /a  Ta  l^^37 

68  L  77  67  68'  78-1 

a  ra  -  (a  /a  )a  T 

-'  (a>;yTa 

68  L  77  '  67  68  78-I  37 


b  =  -a  + 

24  24 


b  =  a  -t- 
33  33 


b  =  a 
34  34 


+  a  ) 


a 

+  a  ) 
^68 


a..) 


a  38^ 
68 


b  =  a  —  a 

55  55  a  51 

11 


pa  -  (a  /a  )a 

l,„  .15  „■  _ Li. _ ±Jk  (a.  -t-  a  ^ 

^a  ^‘^Sl  52^ 

L  41  J 


b  =  a  +  a 


(a  +  a  ) 


55  56  a  '  51  52' 

4  1 


b  =  -a  +  — ^  a 

83  83  a  88 

68 


!a  -  (a  /a  )a  1  a 

_Z2 _ ,,  63^  .6,8  .7.8  (a. - ^  a  ) 

a  -  (a^  /a  Ja  87  a  88' 

77  67  6  8  T'-B-*  6  8 


a  ra  -  {a  /a  )a  I  a 

b  =  a - ^  a  -  -2^ - .  A<.  >6  JL§  (a  -  ^  a  ) 

84  84  a  88  |a  ^  (a  /a  Ja  87 

68  L  77  67  68  78^  68 


b  =  a 

85  85 


a  ra  -  (a  /a  )a  I 

-  -iS  a  -  Lsr,  .LL-.^L  (a  -ha  ) 

a  81  f  2a  J  SI  82 

11  •-  41  -* 


b  =  a  +  a  -  (a  +  a  ) 

86  85  86  a  '  81  82 

4 1 


77  78 

08 
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b 

b  I  ^^b 
a6  b  25 
55 


b 

23  0  24 

34 


b  -T-^b 
86  b  85 
55 


b  -T-aab 

83  D  84 
34 


.  b  b  b  b  . 

[*]  r  =  -2a  a  a  b  b  b  f(b  “vT^^b  )  (b  -Tr^ab  )”(b  ~T.'^'^b  )  (b  -y-5feb  )1. 

11  11  41  68  34  55  77  1.  26  b  25'^  83  b  84  23  b  24  86  b  Bj5  J 

a-  *5«c  f>A.  *a4. 


By  a  similar  diagonalizatlon  process,  ve  can  find  each  of  the  11^^11  ;  e«g*. 


(A. 34) 


"^0,0 

a 

11 

0 

0 

a 

-a 

15 

0 

0 

0 

a 

22 

-a 

23 

-a 

24 

a 

25 

-a  -a 

26  27 

a 

28 

0 

0 

a 

33 

a 

34 

0 

0 

a 

37 

-a 

38 

0 

-a 

41 

0 

0 

a 

45 

a 

45 

0 

0 

0 

-a 

52 

0 

0 

a 

55 

a 

56 

0 

0 

0 

0 

a 

63 

-a 

64 

0 

0 

a 

67 

a 

68 

0 

0 

a 

73 

-a 

▼  4 

0 

0 

a 

.77 

a 

78 

0 

-a 

82 

-a 

83 

a 

84 

a 

85 

a 

86 

-a 

87 

-a 

88 

MRD  DIVISION 

OENERAI.  AMERICAN  TRANSPORTATION  CORPORATION 


A-15 


where 


c  =  a  -  a 

55  55  a  52 

41 

a  a  -  (a  /a  )a  a 

c  =  -a  +^a  -  -  4  5^  4.1,  5.2,  Ua  +,  -AS  a  ) 

26  26  a  22  'a  -  (a  /a  )a  '  ^25  a  22' 

41  55  '  45'  4l'  52  41 

a  a  -  (a  /a  )  a  a 

c  =  a - a  -  (-^6 - -  -*5.  ■  4„J>  .52)  (a - 45  ^  ) 

86  86  a  82  'a  -  (a  /a  )a  '  '  es  a  e* 

41  55  45  41  52  41  ^ 


M=  p  a  a  b  c  b  [p  (b 

0,0  41  68  34  55  77*^*  26' 


b  b 

^  b  )  -  c  (b  -  ^  b  )  ] 

b  84'  86'  23  b  24 

34 


83  ‘J  84 

34 


(A. 35) 


A  =  - 

1 


Thus  after  simplifying  we  find  from  equations  (A. 33)  -  (A. 35)  that 
[c  (t  b  -  b  b  )  -  c  (b  b  -  b  b  )  ] 

OjO  55  .26'  83  34  33  84  86  23  34 _ 33  24  _ 


.55  56l. 


4  33  84  86  55  56  85  23  34  33  24' 


(A. 36) 
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f 

I 

I 


Substituting  Iztto  the  definitions  of  the  b  and  c  coefficients  the 


values  of  the  ve  have: 

m  h  q  n 

b  -  -p  e  *  — «-« - 

*»  *  M  (n  ^+4=) 


m  h  (m  +n  )h  -n  h 


MS  4A  a  I  a*  /  AS  .  —as  .  as  .  ■  m 

™  n  m  2  2_a  2  2'  1  22] 

®  ri  Ch'  -K  ■)■  ■  n  ■)■  -J 

A2ai_g  212J 


-m  h  q  m  n  r 

b  -  -p  e  *  1  -  -g  a  8 -  e 

•*  «  M  (n/+4^) 


„  V,  .  -m  h  (n  -ra  )h  -n  h 

““1  n  -  22-  2  21  2S 

2  1  .  o  ^  .  .  .  -  .  . 


n  (h  - 
»  2'  2 


aaa  y  aa  aa  .  aa  .  a 

k  a  21  22 

i'T' n~CK -VT' 

I'.e  2'  1  s'  J 


n  h  m  h  p  u  (n  ®+f®)-q  n  m  h  -m  h 

^  e  ^  ■  -So-  w; - ^)(e  ‘  ^-e'  ^ 


2S  1'  1 


n  h  -n  h  p  u  (n  m  h  -m  h 

b  -  q  n  (e  ^  i-e  1  1)  -  - -  (e  ^  ^-e  »  M 

26  ^1  1'  q^m^  ' 


b  ap  e 
33  z 


2  r 

h  q  m  n  . 

2  2^  ^2  2  2 _  g 

L 


„  V,  h  (n  +m  )h  -n  h  n  h  2n  h  -n  h  . 

(in^+n2)h^-n^h^  (e  2  z_g^  z  2^  1  =^)re  ^  ^+e  ^  ^  ^ 


-m  h  q  “"m  n 
b  «p  e  ^  ^  .S - 

M  (n  ^+e=) 

z'  z  ' 


n  h  2n  h  -n  h 
e  2  ®  -e  2122 

™  w  -m  h  (n  -m  )h  -n  h  n  n  'dn  h  -n  h 

^n^-ra^jh^-n^h^  (e  *  *-e  *  *  *■  *  ®)  (e  ®  «+e  *  ^  ® 

•"  fTT  iJn  h  -  n  h 

e22_g  21  22 


n  h  mb  -m  h  q  m  n  mb'  -mb 

n  b  -n  b  mb  -mb 

^56  “  ^  ^+®  ^  ^  ^-®  ^ 

n  b  2n  b  -n  b 
®  ®-e  2122) 


TT  2  2 


%[i  - 


1  m  b 
ije  *  1 


-m  b 


M2(n2®+S*)J 
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nh  qmn  mh  m  (n  h  -m  h 

-  «e  *  ‘  e  *  ‘  '  (e  ^  •  *  *) 


V  -  «(' 


n  h  -n  h  M  (n  ^+i^)  m  h  -m  h 
1  1  +  e  1  1)  -  -a— 1- - l(e  1  1  +  e  ^  i] 


ok  - J 


„  „  n  h  -m  h 

°  ^  1  -  e  ^  i) 


ss  1  '  1 


1  _  e  »•  1  x] 

n  h  -m  h 

j  1  1  _  e  -I 


(A. 37) 


Substituting  (A. 37)  Into  (A. 36); 


n  h  -m  h 

^  ^  -  e  1  ^)  X 

I  rp,p,(n /+e=), 


L  o,»,  '  e‘‘>^.e'v*  'J 

FL  (-Jii _ l)e”='’‘l 

[L  =\(n^=+e^)  '  J 


/  ».  .  -m  h  (n  -m  )h  -n  h  n  h  2n  h  -n  h 

(n^-m  )h  -n  h  s  s,^'  e  z'  i  s  ew^  s  e^^  e  i  s  e^^ 


-m  q  n  .  (n  -m  )h  -n  h  /  "  V  e  e  e'"i  e  ew  *e"e 

*^2  ..  /  E.  \\  nTi  tin  ri  -  n  h 

M  (n  *+fi2)'  gEE_e  El  EE 


.[m  (l - ^ - )e 

L  M  (n  M=)'^ 


-m  h  1 
2  1 
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L 


I' 


[(  V  ,  - 1)»  -  f(i  -  - 

lV,(n.*+(»)  •  J  L''  n 


'  )m  e 

An^y  ‘ 

2  2 


_  m  V,  rt  v<  v^  ®  (n  +m  )h  -n  h  n  h  2n  h  -n  h  -Tl 

fp  ».>■.  1.  V. _ /JV".)V"A  (.  »  »-e  °  »  >  »  ■)(;  »  °  »  »  »h] 

r*  -  ">»  ;  'jj 

fr  n  h  -n  h  M  (n  m  h  -m  h  1 

-  M«(e  1  i+e  1  1)-  ■  iq - i(e  ^  i+e  ^  i)J 

r  „  nh  mh  -mh  q^mn  mh  -mhl 

[M,(n^“+«=)(e  ‘  ‘  ‘  ‘)-  -4p^(e  *  *-  *  ‘)J 

nh  -nh  mh  -mhl 
Mjn^®+4®)(e  1  i+e  ^  ^-e  ^  i-e  ^  ^)\ 

f  n  h  q  m  n  m  h  /i  (n  ®+4*)  m  h  -m  h  1 

[{.!>■  .11.  .  X  ^  .  l.)JJ 

fr  _v,  mh(n+ra)h-nh  i 

[.m  vi  fl  2  /  m  Nv,  „  h  “"I  )h  -n  h . n  h.  .  2n  h  -n 

D  ^28  ^2  '^2^2  /-^V"'2^^rVa  (e  g  g-e  ^  ^  ^  (e  «  g+e  «  ^  ‘ 

*^2  ,,  /_  2.42\»*  n  ^  2n  h  -  ri  h 

z 

[m  H  n  -  2™  n  / m  ^ v,  »,  v,  ™  (n  +m  )h  -n  h  n  h  2n  h  -n  h  >, 

P  Va  Va  /.  Va  (e  «  ^-e  =  «  ^  «  «)  (e  «  g+e  «  ^  ^)\] 

*^2  ,,  /„  2j.42\'  n  h  2n  h  -  n  h  l\ 

MgVi^a  ^ 

r  m  >,  «  2_  _  „  .  -rah  (n  -ra  )h  -n  h  -  n  I 

b.'  '  n. (h  !h.)  n.(h.-h  ))J 


.  r,  m  .r.  r.  r,  ^  (n  -m  )h  -tt  h . h  h .  .  2n  h  -n  h  _ 

/^^"2“"‘2^"r"2"2_  (e  g  g-e  2  8^  1  2  2)(e  8  8+e  2  1  2  2)^1 

~  fTTi  2n  h  -  ri  h  }\ 

-22-21.  22  J 


n  h  2n  h  -  n 

e  2  2  _  g  2  1  i 


m  «  2_  _  „  .  -rah  (n  -ra  )h  -n  h 

2  1  <1-  “  "  /  ""i-h  a  2_-'  8  a'  1  2  2 

2  ,,  /_  2.*2\\  n  (h  -h  )  n  (h  -h 

+%  ;  e  ®  ^  -e  2'  i  a 


g  a  a  1  _g  2  1  a' 


(A. 38) 
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similar  expressions  can  be  foimd  for  Aq. 

The  substitution  variables  m^^,  n^,  p^,  can  be  expressed  in  terms  of  the 
original  material  parameters  and  the  variables  ^,s  by 


m  = 
1 


p^(X  +2p  ) 


V  '' " 

/  X  +2jLi 

^  V 


p  =2ju4^+(X  +2u  ) 
1  1  11 


2]Li  +  £S.  (X  +2jLt  )  s‘ 
2  Pi  ^  1  ^i' 


I  (A. 39) 


2p  ^ 
1 


2p  ^ 
2 


The  coefficients  A.  can  then  be  evaluated  and  0  j  ^  ^  ^  ^ 

J  ijOSfOlflSfl 

can  be  immediately  determined  from  equations  (A.15)  -  (A.18). 

A. VIII  Inversion  and  Stress  Solution 


Once  and  ,  are  known,  0.  and©,  are  obtained  through 

1^0  X^X  X  X 


application  of  the  inversion  formulae 


*1-^1  !  ^  *±,o  Jot'-?)®®''  '»=■»? 

Bn  o 
1 


(aAo) 
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